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Non-metricity in the continuum limit of 
randomly-distributed point defects 
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Abstract 

We present a homogenization theorem for isotropically-distributed point 
defects, by considering a sequence of manifolds with increasingly dense 
point defects. The loci of the defects are chosen randomly according to 
a weighted Poisson point process, making it a continuous version of the 
first passage percolation model. We show that the sequence of manifolds 
converges to a smooth Riemannian manifold, while the Levi-Civita con¬ 
nections converge to a non-metric connection on the limit manifold. Thus, 
we obtain rigorously the emergence of a non-metricity tensor, which was 
postulated in the literature to represent continuous distribution of point 
defects. 
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1 Introduction 


The study of defects in solids with imperfections is a longstanding theme in 
materials-science. One of the prototypical crystalline defects are point defects 
(see e.g. lKrd81llKrd901 1. In crystalline materials, point defects may be caused 
by vacancies, interstitials, or impurities. In amorphous materials, point defects 
can be viewed as either a contraction or a dilatation of local equilibrium dis¬ 
tances between adjacent material elements. Assuming that a defect-free body 
is modeled by a smooth Euclidean manifold (see e.g. fKMSlS j). a body con¬ 
taining isotropic (i.e. ball-shaped) point defects is modeled by a subset D c ]R d 
endowed with a Riemannian metric of the form. 


Qr = (p ■ 9 



\x - Xj\ < i for some i = 1,..., m 
\x - xi | > i for every i = 1,..., m, 


( 1 . 1 ) 


where e is the Euclidean metric on ]R‘ f , the points X\,... ,x m are the centers of 
the defects, | ■ | is the Euclidean norm, 1 /R is the radius of a defect, and c is the 
dilatation factor. In this work we focus on defects of vacancy type, hence £ < 1. 
X-, can be thought of as the loci of "missing", or "smaller" atoms in the material; 
"neighboring" atoms thus occupy the vacant location and get closer. 


In mechanics and materials science, a major theme is the modeling of materials 
that contain distributed defects (see e.g. lNol6 7'. fVVan67 il. In continuum mod¬ 
els, bodies with distributed defects are modeled as smooth manifolds, in which 
the singularities are smoothed out (or homogenized) and their density is repre¬ 
sented by an additional geometric field. For example, bodies with distributed 
dislocation-type defects have been modeled since the 1950s as a Riemannian 
manifold endowed with a metrically-consistent, non-symmetric flat connection 
(e.g. 1 Nve53 FBBS551 BS56 1). In this model, the density of the dislocations is 
represented by the torsion-field of the connection. 


A model for distributed point defects has been much more elusive. It has been 
suggested that bodies with distributed point defects could be modeled as Rie¬ 
mannian manifolds with a flat, symmetric, non-metric connection (e.g. llKr68ll 
p. 300-304 and | MR02 |). There is, however, a big difference between the con¬ 
tinuum models of dislocations and point defects. Since the 1950s, there has 
been a clear rationale—even if not a rigorous derivation—relating dislocations 
to torsion. We are not aware of a similar rationale relating point defects to 
non-metricity. In the words of Kroner in his seminal review ( lKro811 p. 304): 
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We are, however, completely aware of the fact that these ... identifications 
[of point defects and non-metricity] have not the same degree of certainty 
as the... identification ... of dislocations and torsion. 


In this paper we present a rigorous analysis of the homogenization of point 
defects. Similarly to the homogenization of edge-dislocations IKM15llKM15bl . 
we obtain a manifold endowed with a non-metric connection as a limit of mani¬ 
folds with distributed defects. Specifically, we consider a sequence of manifolds 
with increasingly-dense point defects. As the density of the defects tends to 
infinity, the sequence of locally-smooth manifolds with singularities converges 
to a smooth Riemannian manifold. The Levi-Civita connection on each man¬ 
ifold in the sequence is, wherever it is defined, the Levi-Civita connection of 
IR"'. When the distribution of the point defects is not uniform, the Levi-Civita 
connection of ]R lf is inconsistent with the limit metric, i.e. has a non-zero non- 
metricity tensor. This is the source of non-metricity in the limit. A surprising 
feature of our result is that the limit metric (and hence the non-metricity tensor) 


is not that expected from volume vs. length considerations, see Section 2.1.3 


In this paper, we investigate isotropic distributions of isotropic point defects. 
That is, the distribution is locally invariant to rotations (isotropy of the distri¬ 
bution), and the defects are ball shaped (isotropy of the defects). A natural 
way to achieve such a distribution is to randomly select the loci of the defects 
using a weighted Poisson point process. The precise model, which is detailed 
in Section 2.2 turns out to be a continuous version of the first passage percola¬ 
tion model on the Euclidean lattice, thus making it an interesting probabilistic 
model on its own; see Subsection |Z2] for more details. 

The structure of this paper is as follows: In Section |Zlj we give a rather informal 
presentation of the main results (without getting into the probabilistic details), 
and discuss their geometric and materials-science/mechanics consequences. 
Section |2T~| is the most relevant for the geometric and material-science-oriented 
reader. In Section 2.2 we describe the probabilistic model for the distribution 
of the point defects and discuss its connections to the probabilistic literature. 
After a list of definitions and notations in Section [23] we state the main results 
in Section 2.4 Since the proofs are rather technical, we provide in Section 3] a 
sketch of the proof of our main results. The detailed proofs are presented in 
Sections 0][8] 
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2 Setting and main results 


2.1 Overview of the results and discussion 

Let d > 2 be the dimension and let E, £ (0,1) be the dilatation factor of the 
point defects. Ignoring momentarily the probabilistic details, our main result 
(Theorem |2.2} is roughly as follows: 

There exists a u, > 0 and a continuous monotonicnlly-decrensingfunction 
p: [0, ii t ) —> (0,1] such that the following holds: Let D c R rf be a compact 
d-dimensional manifold with corners. Let u : D —> [0, iu) he a continuous 
function. Let (D, Qr)r>o he a family of manifolds containing R d ■ J u 
point defects of intensity E, and radius 1/R, randomly distributed in D 
with distribution u. Then (D,Qr) converges (in the Gromov-Hausdorjf 
sense) as R —» oo to the Riemannian manifold (D, (tj o u) • e ). 

This result holds also if E, = 0 (i.e., if point defects correspond to the removal 
of a ball and the identification of its boundary as a single point). The case of 
E, = 0 involves however semi-metrics, hence to simplify the presentation we 
will only consider in this subsection the case E, > 0. 

The next subsections discuss geometric and material-science consequences of 
our main theorem. 


2.1.1 Non-metricity 

The Riemannian (Levi-Civita) connection for each of the manifolds (D, g«) coin¬ 
cides with the Euclidean Levi-Civita connection V of IR‘ ( , whenever it is defined, 
i.e., everywhere except for the boundaries of the defects. Thus, as R —» oo, the 
connection converges (in L°°) to the Euclidean connection. If u is not constant 
(i.e. the point defects are distributed non-uniformly), then the limiting Rieman¬ 
nian metric is not Euclidean, hence the limit connection does not coincide with 
the Levi-Civita connection of the limit metric. In other words, parallel transport 
with respect to the limit connection V is not an isometry in the limit manifold. 

In fact, one can consider the convergence of Riemannian manifolds with con¬ 
nections, 

(D, q r , V) -4 (D, (tj o u) ■ e , V), 

in which case Riemannian manifolds with metric connections converge to a 
Riemannian manifold with a non-metric connection. 

If i] o u is differentiable, then the non-metricity tensor Q; a) of V is given in 
coordinates by 

, , dUn o u) 

Qkij = VkQij = dk&ij - T\ k Qij - r l jk q n = d k (rj o u) <5 i; - = p ^ g /; -. (2.1) 
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where g,-y are the coordinate components of the metric (/] o u) ■ e, and are the 
Christoffel symbols of V, which are identically zero. 

In particular, the non-metricity tensor is diagonal with respect to the metric. 
This is consistent with the model presented in I YG12 1 for bodies with dis¬ 
tributed point defects: Riemannian manifolds with flat, symmetric, non-metric 
connections with a non-metricity tensor diagonal with respect to the metric (in 
IIYG12I such manifolds are called Weyl manifolds). 

We believe that the fact that the off-diagonal components of the non-metricity 
tensor are zero is only a result of our choice of isotropic point defects (i.e. ball¬ 
shaped). For different choices of point defects we expect the emergence of 
non-diagonal non-metricity tensors; see open questions below. 


2.1.2 Curvature 

An immediate corollary of our main result is that any subset of 1R" endowed with 
a conformally-flat Riemannian metric can be obtained as a limit of Euclidean 
manifolds with point defects. 

In particular, the limit manifold can have non-zero curvature, even though the 
point defects do not carry any curvature charge. This is similar to the case 
of dislocations, where it is only the limit connection that is flat—not the limit 
metric (see lKM15b l). 


2.1.3 Length-volume inconsistency 

The Gromov-Hausdorff convergence of (D, Qr) to (D, (ij o u) • e) as R — » oo is 
a convergence of distance functions in metric spaces. Another property that 
converges as R — > oo is the measure Vr induced by the Riemannian metric Qr. 
It weakly converges to the measure p 0OU induced by the Riemannian metric 
(a o u) ■ e on D, where 

o(u) = (e- UK * + £ d (l - e- UK “)) 1/d , 


and Krf = Y(dp~+i) ' s the volume of the d-dimensional unit ball. Indeed, for 
constant distribution u, the defects cover as R —» oo a fraction of (1 - e~ UKd ) of 
the manifold; this remains true locally for non-constant u. In other words, the 
sequence of metric measure spaces (D, dist^, u r ) converges in the measured- 
Gromov-Hausdorff topology (see Definition |2 .1 1 below) to the metric measure 
space (D, d® ou ,p aou ), where dist^ and d^ ou are the intrinsic distance functions 
induced on D by Qr and (/] o u) ■ e, respectively. See part 2 of Theorem 
Corollary |2.3| for details. 

A naive guess would be that ij equals a, since then, both the limit distance 


2.2 


and 
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function and the limit measure are derived from the same Riemannian metric 
(i a is the d -th root of the volume reduction), like the distance functions and 
measures for every finite R. Even though our analysis does not yield an explicit 
formula for 77 , we show that, in fact, i]<a (part 3 of Theorem |2.2j l and in the case 
£ = 0, even 77 < a d (actually 77 < a d is achievable, see a remark in Section 


6.21. 


Hence, the limit metric and the limit measure are inconsistent with each other. 
This can be viewed as another type of non-metricity, not related to the connec¬ 
tion, which, to our knowledge, has not been mentioned in the material-science 
literature so far. 


2.1.4 Open questions 

We conclude this section by raising several natural questions awaiting further 
analysis: 

1. In the present work we assume that point defects are spherically-symmetric 
("isotropic" point defects). A natural question is what is the limit if one 
takes non-isotropic point defects, say ellipsoids. In the non-isotropic 
case, our analysis predicts convergence to some limiting metric space; 
the latter is not conformally-Euclidean, unlike the isotropic case. It is 
not clear, however, whether the limit distance function is induced by a 
Riemannian metric (a plausible alternative would be a Finsler metric). 

If the limit distance function is induced by a Riemannian metric, then we 
expect the resulting non-metricity tensor not to be diagonal with respect to 
the metric as in l |2.1 ) . Either way, whether the limit metric is Riemannian 
or Finsler, this will show that Weyl manifolds (in the sense of [ YG121 1 are 
not the most general model for distributed point defects, as suggested in 
liYGl2l . 

2. A similar question arises if the defects are placed on a grid determinis¬ 
tically (or if the distribution is not isotropic). As the grid spacing tends 
to zero, we conjecture the appearance of a non-Riemannian limit metric, 
even if the entire structure is symmetric (say, cubic defects on a cubic 
grid). 

3. This work focuses on the phenomenological question of describing bodies 
with distributed point defects. Another natural question is how do the 
non-metricity and the length-volume inconsistency manifest in the me¬ 
chanical, or elastic properties of the body. That is, if each manifold (D, g«) 
represents an elastic body with some elastic energy density related to its 
metric, what is the elastic energy functional in the R -4 00 limit? 
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2.2 Probabilistic setting 

Let d >2. We consider the space of locally finite point measures. 


n = 



X; £ for all i > 0 and co{A) < oo 
for all compact A C ]R rf 


( 2 . 2 ) 


with its natural cr-algebra generated by the evaluation maps co i—> co(A), with 

A running over all Borel-measurable sets in !R". 

For x e 1R‘ / and r > 0, let 


B(x, r) = {y e R d : \x - y\ < r) 

be the closed Euclidean ball of radius r around x. 

Let £ £ [0,1). Given co £ Q and R > 0 we denote 

S R (co) = IJ B(x,l/R), 

xesu.pp(co) 

and define a Riemannian (semi-)metric on lR rf 


q r (x; co) = 



x t S R (co), 
X £ Sr(<o), 


(2.3) 


where e is the Euclidean metric on 1R‘ ? . We will often remove co from the notation 
when no confusion occurs. For £ > 0, Qr is a Riemannian metric, and for £ = 0 
it is a semi-metric. Let dish be the distance function induced by g r, that is 

dish(x, y;a>) = inf{lenK(y) : yeT(x, y)), (2.4) 

where T(x, y) denotes all the paths from x to y, and len/dy) is the length of y 
induced by g«. If £ = 0, dish is a semi-distance. 

A note on nomenclature: the term "metric" is commonly used in two different 
contexts—for a Riemannian metric on a smooth manifold and for a distance 
function in a metric space. Since the distinction between the two is at the heart 
of the present work, we will consistently call the first a metric and the second 
a distance. 

Finally, denote by vr the measure on 1R :,; induced by gR, 

vr(A; co) = Leb d (A \ S R (to)) + £ d Leb d (A n S R (co)), (2.5) 

where Leb d is the d-dimensional Lebesgue measure and A c R rf is a Lebesgue 
measurable set. 
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The triple (R d , dish, v r) is a metric measure space if 5 > 0. In order to obtain a 
metric measure space for 5 = 0, we define the equivalence relation, 

x r ~' : y <=> x = y or x, y are in the same connected component of Sr(co). 

_( 2 - 6 ) 

In other words, for every x € supp(cc), we identify all the points in B(x, 1/R). 

For given co e Q and R > 0 the equivalence relation yields a metric measure 
space (Mr, distR, vr), where 


Mr = U d /~. (2.7) 

Denoting by 7 Tr : ]R rf —■» Mr the equivalence class map associated with ~ (and 
with a slight abuse of notation), 

dist R (x, y) := dist R (n^ 1 (x), n~^(y)), Vx, y e Mr. (2.8) 

and 

v r (A) := v R (n^(A)), V Borel set A c Mr. (2.9) 


Note that in | |2.8) both n^ 1 (x) and n~^(y) may contain more than one element. 
However, the distance doesn't depend on the choice of the representatives by 
Definition IZ41 of distR. We denote by M' R the set \x e Mr : Ire)) 1 (x)| = 1}, where 
here and below | ■ | denotes the cardinality of a set. 


Comment: In simple words, given w e Q and R > 0 the metric measure space 
Mr = Mr(o>) is the metric measure space induced from IR" by identifying all 
the points in the closed balls B (x;, 1 /R), where x, are the points in the support 
of co. 

Similarly to the full space, for a path-connected Lebesgue-measurable subset 
D C ]R rf , we denote by d ist^ the intrinsic distance/semi-distance induced by Qr 
on D, that is 

distj^ = inf{lenR(y) : y e T(x, y), y c D), (2.10) 

and denote by the restriction of Vr to D. For £, > 0, this yields a metric 
measure space (D, dist°, v°). For 5 = 0, we denote 

Dr = D/ (2.11) 

obtaining a metric measure space (Dr, distj^, v I ! !), where dist{( and are defined 
(with a slight abuse of notation) as the pullback of dist° and v° by 7Tr, similarly 
to the definitions in | |2.8) and | |2.9) . 

Note: in order to address at the same time the cases 5 = 0 and 5 > 0, we will 
sometimes write Dr and tir even if 5 > 0. In this case Dr = D and 7Tr is the 
identity map (when 5 > 0, all points in D are distinct). 
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Given a function u : R d —> [0, oo) and R > 0, we denote by P U/ r the probability 
measure on (Q, T) under which a; is a Poisson point process with intensity 
R d u(x) ■ Leb d(dx), see HRes871 Chapter 3] for details on Poisson point processes. 
We denote by E U/ r the corresponding expectation. Note that as R —> oo, the 
number of point defects grows like R d whereas the volume in R rf of each point 
defect scales like 1/R d , which is why we expect the measure to converge. 

We will show below (Lemma |4.1) that for £. = 0 and for every density function 
u : R rf —» [0, oo) taking values in a compact set of an interval [0, u*), u » = iu(d) > 
0, the metric space Mr is P„ ^-a.s. simply-connected and locally isometric to 
the Euclidean space, everywhere except for a nowhere dense set. This is also 
trivially true for the case £, > 0, in which (R rf , g r) is locally isometric to the 
Euclidean plane at every continuity point of Qr. 


2.2.1 Discussion on the probabilistic model 

If the function u : R rf —» [0, oo) is spatially invariant (i.e. constant), one can 
view our model as a continuous version of first passage percolation on Z rf , in 
which one associate a random weight with each edge independently according 
to some weight distribution F, see for example lKes86l . Our specific choice of 
points defects is related to the function F(c) = p and F(l) = 1 — p for a certain 
choice of p that depends on it. Note however that one can easily generalize the 
model to general weights by sampling the dilation factor of each ball randomly 
according to some distribution. 

The discrete model of first passage percolation was introduced by Hammersley 
and Welsh 1HW65 1. and was studied extensively since then; see e.g. I Kes86i 
1GK841 lKes93l IBKS03 1 and the references therein. The continuous version of 
this model is based on its continuous counterpart for percolation (known as the 
Boolean model); see l Hal85i. MR96 I. To the best of our knowledge, there is no 
existing work on a continuous version of first passage percolation. 

For a general continuous function u, the discrete counterpart of the model is 
first passage percolation with independent but not identically distributed edge 
weights. 


2.3 Notation and definitions 

In this subsection we list notations and definitions that will be used throughout 
the paper. 


2.3.1 Norms 

• | -1 - The inner-product (Euclidean) norm on R rf . 
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• I • loo - The sup remum norm on IR‘ / . 

• II ' llco - The supremum norm on continuous functions. 


2.3.2 Riemannian metrics 

• e - The Euclidean metric on IR d . 

• Qr = 9 r(-; a>) - The (random) metric on 1R" induced by point defects ac¬ 
cording to (|2.3) (or equivalently (|1.1)). 


2.3.3 Distance functions 

• Given a continuous function p : R rf —» (0, oo), 

^ _ the distance function induced 
p by the Riemannian metric p(x) • e. 

• For a path-connected compact subset D c ]R rf , 

_ the intrinsic distance function induced 
P by the Riemannian metric p(x) ■ e on D. 


( 2 . 12 ) 


(2.13) 


• dists = distp>(-; co) - The distance function induced by Qr on IR 1 ' or Mr (see 
Tft and (fZ8l). 


• dist - a shortened notation for dish used in Sections [4}|6] 

• distj^ = dist^(-; a>) - The intrinsic distance function induced by Qr on a 
path-connected Lebesgue-mesurable set D c IR‘ ? , or on Dr (see |2T0j and 

B ). 

• For two compact sets A,B c ]R‘ ; , we denote by du(A, B) their Hausdorff 
distance with respect to the Euclidean metric on IR d , 


d H (A,B) = inf < e > 0 : 


Vx e A 3y e B, s.t. |x - y\ < e and 
Vx e B 3y e A, s.t. |x - y\< e 


• For compact metric spaces X and Y, we denote by dcu{X, Y) their Gromov- 
Hausdorff distance (see Definition|2.1|below). 
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2.3.4 Measures 












Lebrf - The Lebesgue measure on R‘‘. 

j.Lp - The measure induced on 1R‘ ? by the volume form of the Riemannian 
metric p ■ e. 

/ip - The restriction of p p to D. 

Vr = vr(-; co) - The measure induced by g r on ]R rf or Mr (see < |2.5) andl 2.91). 
v r = v r("' (y ) " The restriction of v R to a subset D c R rf (or to D R ). 


2.3.5 Other notations 

• B(x, r) - The (open) Euclidean ball of radius r > 0 around x € R d . 

• A - The closure of a set A c ]R d . 

• d iam c (A ) = supj|x — y\ : x,y £ A) - The Euclidean diameter of a set 
A c R rf . 

• S rf_1 - The Euclidean sphere {ieK' i : \x\ = 1}. 

• SO(d) - The special orthogonal group in IR rf with respect to the Euclidean 
inner-product. 

• kj = T yp +1) - The volume of the d-dimensional unit ball. 

• By a path in ]R rf we will always mean a continuous parametrized path 

[0,l]->lR rf . 

• len c (y) - The Euclidean length of a path y. 

• len/d}') = len/dv; co) - The length of a path y induced by 9 r. 

• len(y) - a shortened notation for leni(y) used in Sections |4j|6] . 

• l x * y\ ~ The linear segment connecting x and y in IR rf . 

• Lip(/) - The Lipschitz constant of a continuous function / : D —> R for 
some D c R" (with respect to the Euclidean metric). 

• C(D) - The space of continuous functions on D. 
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2.3.6 Measured Gromov-Hausdorff convergence 

'Definition 2.1 In the following (Z, dz) is a compact metric space and (Z, dz, pz) is a 
compact metric measure space. 

1. For a function f : ( X,dx ) —> (Y,dy) between tzvo metric spaces, we define the 
distortion off by 


dis/ = sup \d x (x,y) - d Y (f(x),f(y))\. 


(2.14) 


x,yeX 


A function f : D(f) c X —* Y is called an e-approximation if dis/ < e, 
B x (D(f), e) = X and B Y (Im(f), e) = Y, where B X (A, e) is the e-neiborhood in X 
around a subset A, with a similar definition for neighborhoods of sets in Y. 

2. The Gromov-Hausdorff distance is a distance function between isometry 
classes of compact metric spaces. For the purpose of this paper, it is enough 
to state that dGH((X,d x ),(Y,d Y )) < 4e if there exists an e-approximation f : 
D(f) c X —> Y. For further details see e.g. HPet061 . Chapter 10. 

3. A sequence ( X n , d n , p n ) of compact metric measure spaces converges to a compact 
metric measure space ( X,d,p ) in the measured Gromov-Hausdorff topology 
if there exists a sequence f n : D(/„) c X„ —» X of e n -approximations, with 
e n —> 0, such that thepushforward measures ( f„)#p„ weakly converge to p. That 
is, for every continuous function ® : X —»]R, 



2.3.7 Remark about constants 

Throughout the paper, constants are denoted by C and c. The dependence of 
constants on parameters will be denoted by brackets. For example, C = C(d) 
implies that C only depends on the dimension d. Note that the value of such 
constants may change from one line to the next. Numbered constants C\, Ci,... 
have a fixed value which is determined in their first appearance. 

2.4 Main results 

Our main result is the following: 

‘Theorem 2.2 Let d > 2. There exists a real number u, > 0 depending only on 
d, and there exists a continuous, non-constant, monotonically-decreasing function 
/] : [0, uf) —» (0,1] with 77(0) = 1 that depends only on d and S, e [0,1), such that for 
every compact d-dimensional submanifold D c ]R rf zvith corners and every continuous 
function u : D —> [0, «,): 
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1. Metric convergence: For every e > 0, 


lim 

R—»oo 


u,R 


sup \d nou (x,y) - dist R (7t R (x),7i R (i/))| < e 

x,yeD 


= 1 . 


(2.15) 


2. Measure convergence: For every e > 0, 


R —»oo 

where 


lim P„JV/6W(D), 


f f d t l aou~ f f ° n R dv‘ 
Jd Jd r 


< £ = 1 , 


W(D) = {/ e C(D) : H/IU < 1 and Lip(f) < 1}, 


and 

o(u) = (e~ UKd + £ rf (l - e- UKi j) 1,d . 

3. Length-volume incompatibility: 

T](u) < e~ UKi + £(1 - e~ UKd ), 

and in particular tj(u) < o(u); iff, = 0, then q(u) < o d (u) < o. 

4. The following is of interest when £ = 0: for every e > 0, 

lim P U/R (d H (D, n^(D' R )) < e) = 1, 

R—»oo ' ' 


(2.16) 


(2.17) 


where D' R = {x e Dr : In^x)! = 1}. In other words, n R : D' R —> D is 
asymptotically surjective. 

Theorem |2.2| implies the following corollary: 

Coro [[ary 2.3 

1. Part 1 ofTheorem \2.2\ implies that for every e > 0, 

lim P U , R (d GH ((D r , dist£), (D, d° j) < e) = 1. (2.18) 


Moreover, for every sequence ( a) n ,R n ,e n ), where co n e Q, R n > 0 and e n —> 0 
such that the events in pd5)-||2T7) hold with e n and R n , the sequence of 
metric measure spaces ( D Rn , dist^,!'® )(or„) converges in themeasured-Gromov- 
Hausdorjf topology to the metric measure space (D, d l tRJU , pf ou ). 


2. Part 3 in Theorem 


2.2 


implies that while the measures p R and the distance 
functions dist^ on D Rji are consistent, the limit measure p aou is strictly larger 
than the measure p, ]OU induced by the metric (/] o u) ■ e associated with the limit 


distance d 


TjOU‘ 
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Comments: 


1. The constant u t is the percolation threshold; for u > iu, S] contains almost 
surely a unique infinite connected component (see also a remark at the 


end of Section [4~2| . 


2. The assumption of u being continuous can be relaxed; for example, the 
theorem holds if the d-dimensional Hausdorff measure of the set of dis¬ 
continuous points of u is zero. 

3. The bound < o d when c = 0 is not tight. With a more complicated 
argument it is possible to show r\ < o d . See a remark at the end of 
Section 16. 2 1 

4. Throughout this paper, we assume that "manifolds with corners" do not 
include cusps. This assumption can be presumably relaxed. 

3 Sketch of the proof 

In the remaining sections of this paper, we prove Theorem |2.2| and Corollary |2.3| 
For the sake of brevity, we assume that 5 = 0, which is the most difficult case; 
the proofs for 5 > 0 are similar, and in certain parts, much simpler. Since the 
proof is long and technical, this section describe its main stages. 

In Sections [4jj7] we consider a uniform distribution of point defects, that is, the 
case where u is a constant function. In Section [8| we generalize the results to 
arbitrary continuous functions u : D —> [0, u,). 

Section]!] For a constant u, a uniform rescaling of ]R rf enables the introduction 
of a natural coupling between the probability measures (P U/ r)r>o using a single 
measure P„ that is used in the rest of the proof for the uniform case. 

We show the existence of a sub-critical regime, i.e., a constant u, = «,(d) > 0, 
such that for u € [0, it,), the process P„-a.s. does not percolate. In this regime, 
Mr is "nice"; it is P„-a.s. a simply-connected metric space, which is locally 
isometric to the Euclidean space but for a nowhere-dense set. From this point 
onward we work only in the subcritical regime u < iu. 

We then use the subadditive ergodic theorem to obtain the existence of a limit 
distance function in Mr. We show that for every u e [0, iu) there exists an 
7](«) e [0,1) such that 



Moreover, this limit is uniform in every compact K c 1R‘ / . This establishes the 
P„-a.s. Gromov-Hausdorff convergence of (Kr, dish) to (K, d, 
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Note, however, that this is not quite the metric convergence we want in 1 2.15) for 
a uniform u, since we want to prove that (Kr, dist°) converges to ( K , d^ u) ) (i.e. the 


intrinsic distances converge and not only the induced distances). Moreover, we 
need to prove that rj(n) > 0 for u e (0, u t ). This is done in the next sections. 


Section [5j In this section we prove large deviation results for the distance 
function in Mr. An immediate corollary is that r\(u) > 0 for u e (0, u t ). The key 
idea here is to exploit the independence structure of the Poisson point process, 
manifested in the BK inequality, in order to show that distances do not deviate 
significantly from their expected value. 


Section [6] In this section we prove several results regarding the geometry of 
geodesics in (Mr, dist R ), and properties of the function i] that controls the limit 
metric. 

By using the concentration results of Section[5j we prove that geodesics in Mr 
are, with high probability, very close to straight lines between their endpoints. 

Using a coupling between the probability measures (P u )u>o, we prove that r\(u) is 
a continuous, monotonically-nondecreasing function, and give an upper bound 
on the value of //(«) that implies that r] < a (which proves part 3 of Theorem |2.2) . 

Finally, we use the ergodicity of the model to prove that Vr([0, l] rf ) converges 
P„-a.s. to u„(„)([0, l] rf ) (which is a step towards proving the measure convergence 
in part 2 of Theorem|2.2[l. 


Section[7j In this section we prove Parts 1 and 4 of Theorem |2.2| for uniform 
distributions on a convex, compact d-dimensional manifold with corners D C 
]R rf . The idea behind the proof of Part 1 is as follows: from Section [ 4 ] we 
know that P„-a.s. (Dr, dist R ) Gromov-Hausdorff converges to (D, d,j(„)). Since 
D is convex and u is constant, d,,( u ) = dj^ on D, so we only need to replace 

(Dr, dist K ) with (Dr, dist R ). We do so by showing that the identity mapping 
between (Dr, dish) and (Dr, dist R ) has P„-a.s. vanishing distortion. Here we 
use the result from Section [6] that geodesics in Mr are very close to straight 
lines, which implies that distR-geodesics between points tend to remain within 
D, i.e. they are also distR-geodesics with high probability. 


Section |8j In this section, we conclude the proof of Theorem |2.2| and Corol- 
lary |2.3| for a general continuous distribution u of defects over a rf-dimensional 
manifold with corners D C ]R d . 

The idea is the following: we partition D into small cubes, such that u is 
approximately constant in each cube. We show that the results of Sections[6]and 
[7]apply approximately to each of the cubes, in the sense that for every cube □, 
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when R is large, the distance between the metric measure spaces (Dr, dist^, vy) 
and (□, d° oti , is bounded with high probability by the variation of u in □. 

Finally, we glue the cubes together and obtain results for the whole manifold 
D. For measure convergence, the gluing is straightforward. For metric conver¬ 
gence, we use the control on the distortion in each cube and a bound on the 
number of cubes each geodesic crosses to control the Gromov-Hausdorff dis¬ 
tance between (Dr, dist^) and (D, d^) (in a similar way as in IKM151lKM15bi l. 

4 Uniform distribution of point defects 

In this section, as well as in the three to follow, we study the simplest version 
of the model: we assume that the function u : 1R‘ ? —» (0, oo) is constant. In this 
section, we prove the existence of a subcritical regime (Lemma |4.1| , and obtain 
our first main result regarding distances in the manifolds Mr (Theorem |4.2| . 
Many of the results obtained in these sections are adaptations to the continuous 
setting of the results obtained in I 11VV65 Kcs86 l for the discrete case. 

4.1 Rescaling for uniform distributions 

The parameter R > 0 is a scaling factor that affects both the density of the point 
defects and their magnitude, or volume. In particular, to every R corresponds 
a different probability measure P„,r. For spatially invariant u, the intensity 
measures are translationally-invariant. By a uniform rescaling of space, we 
may obtain a probabilistic model that does not depend on R, and in particular 
allows us to construct a natural coupling of the measures P,/,r. 

This is done as follows: construct a Poisson point process on JR‘ ? with intensity 
u ■ Leb d(dx), i.e., set R = 1. We denote the probability measure by P„ = P„j , and 
the corresponding expectation by E„ = E U/ i. Two points in All are identified, 
x a ~ y, if they are in the same ball of radius R = 1 centered at a point in the 
support of a>. We denote by dist = dish the corresponding distance function 
on R d , which can also be written as 



(4.1) 


Given R > 0, let T R : 1 R rf —> ]R rf be defined by T R (x) = Rx. Then, (Ti/ R )#co is 
a Poisson point process with intensity R d u ■ Leb d(dx), namely P u ,r = (Ti/r)#P„. 
Here, (Ti/r)# is the push forward of measures by Ti/r, where the push forward 
of P„ is the one induced by the push forward of a>. That is 



for 
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and. 


(Ti /r )#P u (A) = P u ({w : (T\/ R )#co e A}) 

Similarly, we identify Mr with M\ via the scaling T\/ R . The distance function 
dish on \R d is derived from the distance function dist = disti by the following 
relation, 

dist R (x, y;co) = — dist (m (R n~ R (x)), nfR n~ R (xy))-, or), Vx, y e Mr. (4.2) 

Note that there is a slight abuse of notation here: strictly speaking, dist R (x, y; a>) 
as defined in ( |4.2[ ) coincides with dist R (x, y; (Ti/r)#o>), as defined in l |2.4) . How¬ 
ever, the distribution of ( |4.2| l with respect to P„ is the same as the distribution 
of \2A) with respect to P„, R . In Sections |4j|7j where the above coupling is used, 
we use dist R (x, y; co) in the sense of l |4.2) , so that P„ can be used for all values of 
R. 

4.2 The sub-critical regime 

We start our analysis by proving the existence of a subcritical regime: 


Lemma 4.1 For co e diet S(co) = Si(co) = U^esupp^) !)• Then 

(1) Each of the connected components ofS(co) is P u -almost surely closed. 

(2) There exists a constant u, > 0, depending only on d, such that for all u e (0, u,) 

p / The set S doesn't contain an 
“ I infinite connected component 

and for every u > u, 

p / The set S contains a unique 
" 1 infinite connected component 

In particular, recalling that E, = 0, whenever u e (0, iu) the simply connected 
metric space Mr is P u -a.s. locally isometric to the Euclidean space, up to a 
nowhere dense set. 

Comment: If £ > 0, Mr is almost everywhere locally isometric to a Euclidean 
space, with a scaling constant depending on whether the point is the interior 
of S or in the complement of S. 

Comment: The subset M' R of Mr that is locally isometric to the Euclidean 
space can be identified with IR"' \ S(a>), not only as sets, but also as Riemannian 
manifolds. However, they are not globally isometric. Note that M' R is P„-a.s. not 
connected. 
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Proof: The fact that each connected component of S(a>) is P„-almost surely 
closed follows from the fact that with P„-probability one, every bounded set 
contains only finitely many points of supp(oj). This implies that P„-almost 
surely, the point process does not have accumulation points, hence the comple¬ 
ment of S(co) is P„-almost surely open. 

The existence of u, > 0 such that S(co) doesn't contain P„-a.s. infinite clusters 
for every u e (0,1/,) and contains P u -a.s. an infinite cluster for u > u, is the 
content of [ Hal85l ; see also Theorem 3.3 in I 1MR961 . The uniqueness of the 
infinite cluster for u > u, can be found in IMR96 I Theorem 3.6. 

Since Mr is obtained from IR‘ / by a similarity transformation and an identifica¬ 
tion of points in simply-connected subsets, Mr is simply connected. As proved 
above, IR‘ ( \ S(a>), which is identical to M' R up to a similarity transformation, is 
open. It follows that M' R is locally isometric to Euclidean space. 

Finally, we need to show that for u e (0, u,) the set Mr\M' r is nowhere dense. 
This follows from the fact that any compact subset of ]R rf contains P u -a.s. only 
finitely many points in supp(&;), hence every compact subset of Mr contains 
only finitely many points in Mr\AV r . ■ 

Comment: It can be shown that for u > u, and d = 2, in every box of sufficiently 
large radius R, the distance of every point from the unique infinite component 
is at most O(logR). As a result, the limiting distancs dist R between every pair 
of points is zero. 

4.3 Distances in Mr in the uniform case 

The main result of this section proves the existence of a limit distance function 
in Mr. The precise statement is as follows: 

‘Theorem 4.2 For every n e [0, u t ) there exists an i](u) e [0,1] such that 

P„ (lim dist R (7T R (x), n R (y)) = /](«) ■ |x - y\, Vx, y e = 1. (4.5) 

\R—*oo ) 

The limits also exist in the L 1 ( Q, P„) sense: for every x,y e ]R rf , 

lim E„ [|dist R (7T R (x), 7i R (y)) - rj(u) ■ |x - y||l = 0. (4.6) 

R —>oo Ll 

Furthermore, the convergence of the distance function is uniform in every compact 
Kc~R d (and in particular in S rf_1 ; this particular case is used belozv), 

Ve > 0 3Ro such that VR > Rq, 

SU P V ye R |dist R (7t R (x), n R (y )) - p(u) ■ |x - y|| < e 

and 

lim sup E„ ||dist R (7i R (x),7i R (y)) - ?](») ■ |x - y||l = 0. (4.8) 

R ^°°x,yeK L J 
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In Sections |5jj6]below we prove certain properties of r\(u), and in particular that 
Tj(u) > 0 for every u e [0, m). 

Equation ( |4.5[ ) and similar equations hereafter should be interpreted as follows: 

P u ^jo; e Q : Km dish(7tR(x), Hr(j/); cj) = q(u) ■ |x - y\, Vx, y e ]R rf | j = 1, 

where dish is understood as in ( |4.2[ ). 

An immediate corollary of Theorem |4.2| is: 

CoroUary 4.3 For every u e [0, iu) and compact K c ]R rf , the sequence of metric spaces 
( Kr, dish) defined by l |2dl] P H -«.s. Gromov-Hausdorjf converges to (K,d n ( u )), where 
i](u) should be considered as a constant function on ]R lf . 


Note that this convergence is with respect to the induced distances dish on 
Kr and not with respect to the intrinsic metric dist^ on Kr, defined in 1 2.10) . 
Proving the convergence of the intrinsic metric is more involved and requires 
more assumptions on K. This is done in Section [7] 


4.4 Proof of Theorem 14.21 


We start by reformulating Theorem |4.2| using the relation between dist and dish 
and the definition of Ur\ 


‘Theorem 4.4 ("Rephrasing of‘Theorem 4.2) For every u e [0 ,u t ) there exists r/(u) e 
[0,1] such that 


P„ I lim dlSt ^ /R ‘^ = i ](u) ■ \x - y\, Vx, y e R d ] = 1. 

\R->oo K 


(4.9) 


The limits also exist in L 1 (Q, P u ):/or every x, y e ]R rf , 


lim E„ 

R—>oo 


dist(Px, Ry ) 
R 


-'](») • \x - y\ 


= o. 


Furthermore, the convergence of the dist(Px, Ry)/R is uniform over x,y in every 
compact K c R d (and in particular in S' 1-1 ; this particular case is used below). 


The proof of Theorem 4.4 is separated into several parts and starts with the 
observation that the system (Q, 'T, P„) is ergodic with respect to translations. 


Lemma 4.5 For x e ]R rf define t x : Q —» Q by t x (E,>o = L;>o b x ,-x- Then, for 

every x e ]R rf \ {0} the quartet (Q, (F, P u , v x ), defines a translation invariant ergodic 
system. 
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Proof: See for example 1 MR96 I Proposition 2.6. 

Next, we prove the existence of the limit in (|4.9[) for x = 0. 


Lemma 4.6 Let n e [0, u t ). For every y e R rf the limit 

dist(0, Ry) 
pu(y) = lim - - - 

R—>oo 1\ 


exists P u -a.s. and in L 1 2 3 4 (Q, P,,). 


(4.10) 


Proof: If u = 0, then a>(]R rf ) = 0 P„-a.s., i.e., dist(x, y) = \x - y\ with P„-probability 
one, which implies that 


lim 

R—>oo 


dist(0, Ry) 
R 


= I yl. 


P, r a.s. 


Thus the result holds with po(y) = \y\. 

We turn to the case u > 0. For y = 0 the statement is trivial, so fix y e lR rf \{0} 
and define for 0 < m < n 

Y m ,„ = dist(my,«y). (4-11) 

The triangle inequality for dist implies that 

Yo,n < Yo,m + Y m ,n, V0 < m < n. (4.12) 


Thus we are in a good position to use Kingman's subadditive ergodic theorem 
HKin73l . More specifically we will exploit Liggett's version ] Lig85| , which states 
that if (X m „)o< OT<K are nonnegative random variables such that 


(1) Xo,„ < Xo,m + X m , n for all 0 < m < n, 

(2) {X n i.,( n +i)jt : n > 1) is stationary and ergodic for each k> 1, 

(3) the law of {X m m+ j- : k > 1) is independent of m > 1, and 

(4) X 0 ,i has finite expectation, £[Xo,i] < oo, 


then the limit lim,,-,,*, Xo,„/« exists almost surely and in L 1 and almost surely 
equals 


inf 

n> 0 


E[X 0 , n 


,. E[X 0 J 

= lim - < oo. 


jj n—>oo n 

Taking X m n = Y„ h „, (1) is given by l |4.12| . For (2) and (3) note that 
dist(n +z,v + z;a>) = dist(«, v; t z co), V»,z’,z e ]R rf , 


(4.13) 


and therefore Y„ hn = Y o rH - m ° t™. Condition (2) and (3) then hold by the transla¬ 
tion invariance and the ergodicity of the law P„ under the shift i v , see Lemma 


4.5 As for (4), it follows immediately from the fact that Yo,i = dist(0, y) < \y\. 
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Thus, the family {Y mA }o< m<n satisfies conditions (l)-(4) and 

lim n~ l dist(0,uy), 

n—>oo 


which we denote by p„(y), converges P„-a.s. and in L 1 . We then prove that 

dist(0, Ry) 


Pu(y) = lim 

R—>oo 


R 


P„-a.s. 


(4.14) 


Indeed, p„(y) = lim^K, dlst( | 0 ^j R ~ ly) and since dist(|_RJy,-Ry) < ||_RJy - Ryl < Iy|, it 
follows that 

dist(0,Ry) dist(0, |RJy) ^ 2|y| 

R [Ri “ IT' 


This completes the proof. 

■ 

Lemma 4.7 The function p u in Lemma ^L6\satisfies 


Pu(ay) = ap u (y), Vtv e (0, oo), y e R d , 

(4.15) 

and 


Pu(y + z) < p u (y) + p„(z), Vy,z e ]R rf . 

(4.16) 


Proof: The positive homogeneity follows from the existence of the limit since 

dist(0, Ray) dist(0,Ray) 

Pu(ay) = lim ---- = lim a-—-= ap u (y). 

R—> oo i\ R—>oo l\Ct 

For the triangle inequality note that by the translation invariance of P„ 

E„[dist(0,R(y + z))] < E„[dist(0,Ry)] + E„[dist(Ry,R(y + z))] 

= E„[dist(0,Ry)] + E„[dist(0,Rz)]. 


(4.17) 


Diving both sides by R, taking the limit R —» oo and using the L 1 convergence 
of dist(0, Ry) /R gives the required inequality. ■ 


Lemma 4.8 The function p u in Lemma 4.6 is invariant under the action ofSO(d). 

Proof: Since P„ is invariant under the action of SO(d) it follows that for every 
two points yi,y 2 e with the property that there exists 'R e SO(d) such that 
y 2 = Kyi, 

dist(0, yf) = dist(0, yi). 


Thus, p u inherits the symmetries of SO(d) as the L 1 (Q, P„) limit Lemma 4.6 


Coronary 4.9 For u e [0, u t ), the function p u in Lemma 4.6 is of the form 

Pu(y) = T](u) ■ |y| 

for some rj(u) e [0,1]. 


(4.18) 
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Proof: This follows from the positive homogeneity, the sub-additivity and the 
isotropy of p u proved in Lcmmas |4.7| and |4.8| ■ 

Next, we consider general x, ye ]R' f : 


Lemma 4.10 Fix n e [0, »*). Then, for every x,y e R‘ f 


lim 

R —>oo 


dist(_Rx, Ry) 
R 


='](«) • I* - y\ 


converges P u -a.s. and in L 1 . 


Proof: For every x,y e\ R rf and aieQ the relation dist(x, y; co) = dist(0, y - x, t x co) 
holds and therefore 


dist(_Rx, Ry;co) dist(0, R(y - x); r x co) 

R = R ' 

which implies that the limit R —» oo exists P u -a.s., is in L 1 (Q, P„) and equals 
r](u) ■ \x - y\. ■ 


Proof of Theorem 4.4 By Lemmas 4.6- |4.10 dist(Rx, Ry)/R converges as R —» oo 
both P„-a.s. and in L 1 for every fixed pair of points x,y e R rf . It remains to verify 
that the limit exists P„-a.s. simultaneously for all pairs of points x,y E lR lf . 

To this end, let e e (0,2 n) and let (v,)^ =1 with N = \c{d)l e d ~\ be a set of points on 
the unit sphere S rf_1 that form an e/2-net for S rf_1 (c(d) is a const ant th at depends 
only on d). Since the set t is finite, it follows from Lemma ‘ 


4.10 


that 


disi(Rvi, Rvj) 

lim---— = r\(u) ■ |Vi - v;\, VI < i,j < N = 1. 

R—>oo i\ 


Given x, y e 1R‘ # \{0}, there exist 1 < i, j < N such that 

\x - Vj\ < e and |y - Vj\ < e, 

where x = x/\x\ and y = y/|y|. By the triangle inequality. 


(4.19) 


(4.20) 


dist(Rx, Ry) dist(R|x|z>„ R\y\vj) 


R 


R 


< dist(Rx, R|xh,) dist(Ry,R|y|z;/) 


R 

^ [Rx - R\x\vi\ 

~ R + 

<£(l*l + \y\). 


R 

I Ry - %M 
R 


(4.21) 


and 


| n(u) ■ \x - y | - ?](») ■ ||x|i7,' - \y\vj\\ < e{\x\ + |y|). 


(4.22) 
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Combining both estimates, it follows from (4.191 that 


P„ |lim sup 


R —»oo 


dist(iCr, Ry) 


R 


- '?(«) ■ \x - y\ 


< 2e(|x| + \y\), Vx, y e 1R \ {0} = 1 


The cases x = 0, y = 0 can be included by a similar argument, except that no 
approximation for 0 is needed. Since this holds for every e > 0, the limit as 
R —> oo exists P u -a.s. 

To justify the uniformity over compact sets, note that there are only finitely 
many v/s for a fixed e > 0 and that for compact sets we have a uniform bound 
on the Euclidean norm of both x and y. The same applies for the L 1 (Q, P„) 
convergence. ■ 


5 Large deviation results 

In this section we prove large deviation results for the distance in Mr. As an 
immediate corollary we obtain that i](u) > 0 for every u £ [0, u,). 

The main result of this section is the following. 

‘ Theorem 5.1 For every u £ [0, u t ) and every £ > 0 there exists a positive constant c\, 
depending only on d, u and e, such that for large enough R, 

P„ (dx e U d : |xj = R, \ dist(0,x) - p(n)R\ > eR) < e~ ClR . (5.1) 

The proof of Theorem |5 .1 1 follows the ideas developed for the discrete case by 
Kesten [Kes86 |. 

In addition, we will need a large deviation result for the existence of very 
long geodesics in the Euclidean sense. In order to state it we need another 
definition: Geodesics in Mr can be identified with geodesics in IR^ with respect 
to the semi-distance function dist. Such geodesics are highly degenerate, as 
there is nothing that limits their behavior in S(a>). For x, y £ 1R‘ / we denote by 
To(x, y) the set of geodesics between x and y with respect to dist that minimize 
the Euclidean distance inside »S(aj). We will call such paths true geodesics (see 

‘Theorem 5.2 For every u £ [0 ,iu) there exist positive constants C 2 ,C 3 and a » 1, 
depending only on d and u, such that 

(1) For every R > 0, 

/ 3 a true geodesic path starting at 0 \ R 

“y such that len e (y) > R and len(y) < ^ len e (y) )~ 2 ' ' 
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Figure 1: The path on the left is a geodesic in Mr, where the gray regions 
correspond to Sr(o>). Note that the part in Sr has zero length. The path on the 
right is the corresponding true geodesic. 


(2) For every x,y el R rf and for every R > a\x - y\ 

P„ ( 3y e Y 0 (x, y) : len c (y) > R) < C 2 e~ C3R . (5.3) 

As a corollary we obtain: 

CoroCCary 5.3 r\{n) is strictly positive for every u e [0, u »). 


Proof: Let Ar be the event in l |5.2[ . By Part (1) of Theorem 5.2 

V u (Ar) < C 2 e- C3R , VR > 0. 

Since the Euclidean length of a true geodesic connecting 0 and Re i is at least R, 
for co e A c r , 

1 

dist(0,.Rei) > — R, 


Thus, for large enough R 


E u [dist(0,Rei)] > E„[dist(0,Kei) • 1m] > -R ■ P U (A C R ) > ^-R, 

R a 2 a 


hence 


= lim > i_ > 0. 

r —*oo R 2a 
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Comment: A conjecture in percolation theory is that 

lim i](u) = 0. 

For d = 2, an adaptation of IIGri991 Lemma 11.12] to our setting will show that 
this conjecture holds. In higher dimension this a famous open problem. 

The proofs of Theorem [5T] and Theorem |5.2| arc quite technical. The main idea 
behind the proof is to exploit the independence structure of the Poisson point 
process, manifested in the BK inequality, in order to show that geodesics whose 
lengths deviate from the expected distance in the sense of ( |5.1) , contain suffi¬ 
ciently many disjoint sub-paths (which are roughly independent) whose total 
length deviates significantly from its expected value. Such an event is highly 
unlikely due to large deviations results for independent random variables. 

Since the proofs in this section are technical and since one can use the theorems 
as "black boxes" in the rest of the paper, the reader might wish to skip the rest 
of this section in a first reading. 


5.1 The BI< inequality 


In this subsection we state the continuous version of the well-known BK in¬ 
equality for product measures. For this, we need some additional definitions. 

There is a natural partial ordering on Q, which we denote by <, under which 
co < co' if and only if supp(o>) C supp(aj'). Using it one can define increasing and 
decreasing events in T. An event A 6 T is said to be increasing (respectively 
decreasing) if for every co < co', co e A implies co' e A (i.e., A is closed under 
increasing support). 

For any bounded Borel set Y c ]R J , define the set 

coy = supp(o;) fl Y, 
and for co e Q and Y as above let 

L(co, Y) = {co' e Q : supp(aj) n Y c supp(a/) n Y]. 


In words, the event L(co, Y) contains all configurations that inside Y are larger 
than co. We say that an event A is an increasing event on Y if co e A implies that 
L(co,Y) c A. 


:Definition 5.4 Let A and B be tzvo increasing events on a bounded Borel set Y. Then 


AoB 


there are disjoint sets V, W c ]R rf such that 
< co e Q : V and W are finite unions of rational , 

cubes and L(co, V) c A, L(co, W) c B 


where by a rational cube we mean an open d-dimensional cube with rational coordinates. 
When AoB occurs, we say that A and B occur disjointly. 


26 





Esample,: To illustrate this definition, let d = 2, let Y = [0, l] 2 , and let A be the 
event that there exists a path in Y connecting the left boundary of Y to its right 
boundary, whose length is less than x. Clearly, A is an increasing event, as 
increasing the support of a> can only shorten paths. For a> e A and V C Y, it is 
generally not true that L(a>, V ) c A; L(co, V) C A only if there exists a path in V, 
such that its length is not more than x — a, where a is the sum of the Euclidean 
distances of the path's end points from the left and right boundaries of Y. The 
event AoA occurs if there exist two disjoint paths connecting the left boundary 
of Y to its right boundary, whose length is at most x. a a A 

‘Theorem 5.5 (‘BTfinecfuaCity) Suppose Y is a bounded Borel set in ]R rf and A,B are 
tzvo increasing events on Y. Then for every u > 0 

P „(A o B) < P U (A)P U (B). 

A proof of this inequality in a more general setting can be found in [ MR96 I. 
Theorem 2.3. 


5.2 Key proposition 


The proofs of Theorem |5d] and Theorem |5.2| exploit the invariance of the law of 
P,, under translations and rotations. This implies that we only need to take care 
of large deviation results for the distance between the origin and points of the 
form Re where e\ = (1,0,...,0). We therefore restrict ourselves to the above 
case and define the following: 

For r e IR let 

H r = {xe ]R rf : {x,ef) = r). 

Given r < s and a path y we write 

H r < y < H s 


if all points of y, except possibly its endpoints, lie strictly between the hyper¬ 
planes H r and H s . 

The following variants on the distance will stand in the core of the proofs. For 
N,M> 0 define the random variables. 


Sm,n = inf 



y is a path from {0} X [0,N] d 1 1 
to Hm such that Hq < y < Hm ) 


and 


Sm,n = inf • 


len(y) : 


y is a path from {0} x [0,N] rf 1 to H M such 
that with the exception of its endpoints, 
y c (0,M) x [-4M,4M] rf_1 


(5.4) 


(5.5) 
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It follows from the definitions of Sm,n and Sm,n that 


Sm,N 5 Sm,N- 


(5.6) 


The main estimate used in the proofs of both theorems is stated in the next 
proposition: 


imposition 5.6 Let (x q (M, N)j ^ and (x q {M, /V) j ^ be sequences of independent 

random variables having the same distribution as Sm,n an d 'sm,n- F° r every tzvo 
integers 1 < N < M < R/2 and every real number x>0, 


P„ (dist(0,.Rei) < x) < 


< 



\ 

< X 

\ 

< X 

/ 


(5.7) 


Proof: Since Sm,n 5 ~Sm,n, the second inequality is immediate. Assume y : 
[0,1] —» IR !? is a simple path from the origin to Hr. We choose a sequence of 
points (xq,Xi,...,xq) along y as follows: Define x 0 = y(0) = 0. Assume that 
{xq,X\ ...,xf) have already been chosen such that X; = y(f,;) with 0 = to < ti < 
... <tq. Then we define 

tq+i = minjf e (t q , 1] : | y(t) - x q \oo =M + N}, x q+1 = y(t q + 1 ), (5.8) 


provided such time t exists. If no such t exists, i.e., \y(t) - x, ; | K , < M + N for all 
t e (t qr 1], then we set Q = q and stop the process. 

It follows from the definition of the points x q that 


Q> 


R 

M + N 


- 1 . 


(5.9) 


Indeed, the distance between (x a , ef) and (x a+ \, ef) is at most M+N and (xn, ef) > 
R - (M + N). 

Next, we analyze the path within the time interval [t q/ t q+ i] for 0 < q < Q — 1. 
Since, by definition, \x q - x^+iU = M + N, it follows that there exist j = j(q) e 
{1,2,.. .,d) and a = o(q) e {-1,1} such that 


(x q+ 1 - x q ,e m ) = o(q)(M + N). 


We introduce the hyperplanes, 

H\ = [x e Ml 1 1 (x, ef) = r\, 
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j = !,•••, d. 


reK, 



so that, in particular, Hj = H r . The path yip is strictly restricted between 
the hyperplanes 


m = h 


/(<?) 


and 


H'(q) = H 


M 

h?/e/(,)>+a(?)(M+N) 


Since the distance between these hyperplanes is M + N, we can find two hy¬ 
perplanes H"(q) and H"'(q) at a distance M one of the other, that are strictly 
between H(q) and H'(q) (see Figure[2j. Further, we denote by [t"(q),t'"(q)] the 
sub-interval of [t q , f (;+1 ] in which the path y contains a unique crossing between 
the two hyperplanes H"(q) and H"'(q ). 

Specifically, let m j( q )(q) be the integer m such that H^ is the hyperplane bounded 
between H(q ) and H'(q) that is the closest to H(q). Explicitly, 




(Xq, £j(q)) 

N 


We then define 


=hS 


and 


2 ^ 0 ?) + !)• 


H'"(q) = H K £ , 

v i' mN+o(q)M r 


and let [t"{q), t"'(q)] C [t q , t i;+1 ] be the subinterval defined by 


t”(q) = maxjf e [t q ,t q +i) ■ y(t)eH"(q )}, 


t'"(q) = min {f e (t”(q), t q+1 ) : y(t) e H"'(q)}. 

Finally, denote by y q the path y restricted to the time interval [t”(q), t'"(q)]. 

We are interested in a bound on paths starting within a (d — l)-dimensional 
box of side length N within the hyperplane Ho, see {53J. We denote by m(q) = 
(mi (q), miiq), ■.., ntd(q )) the unique point in Z rf satisfying 


y(t"(q)) e A(q) = H"(q) n (m(q)N + [0,N) rf ) . 

where A (q) is a (d- l)-dimensional box of side length N within the hyperplane 
H"(q). 

Exploiting all the above definitions, the path segments y q satisfy the following 
properties: 

• The images of y„ in R rf are pairwise disjoint. 

• y q connects A (q) c H"{q) to H'"(q). 

• The path y q lies strictly between H"(q ) and H'"(q), except for its endpoints. 
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Figure 2: Illustration of the times t cl ,t q+ i,t"(q),t"'(q) / the points 

x q , x q+lr y(t"(q)), y(t"'(q)), the hyperplanes H(q), H'(q), H"(q), H"'(q), the box A(^) 
with its corner Nm(q) and the path y q crossing from A (q) to the H"'(q). 
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• The path y q is contained in the box 


B(q) = Nm(q) + [-4M,4M] /W_1 X [0 ,a(q)M] X [-4M, 4M] d_;(?) , 
where we used the fact that N + M < 2M. 

• The total length of the paths y q satisfies 


Q-i 

Yj len 0b) < len(y). 

ij=0 


Fix x > 0, fix Q > - 1 and fix (j{q),a(q),m{q)), q = 0,...,Q - 1. We 

denote by A(j, a, m, x ) the event that there exists a piecewise-linear simple path y 
containing disjoint segments y q C B(q) crossing the box B(q) in the j(q) direction 

from A (q) C H"(q) to H"'(q), with le n (lh) - x - 

Since every true geodesic is piecewise-linear, we conclude that the event 

{dist(0,.Rei) < x } 

occurs only if there exists a Q > - 1 and there exist (j(q),o(q),m(q)), q = 

0,..., Q — 1 such that A(j, a, m, x) occurs, therefore, 

P„(dist(0,Rei) < x) < E E P« {A(j,o,m,x)). (5.10) 

The event A(j, a, in, x) is characterized by the existence of a path whose proper¬ 
ties are specified over disjoint segments. Denote by A q (j, a, m,x) the event that 
there exists a piecewise-linear simple path y containing a segment y q C B(q) 
crossing the box B(q) in the j(q) direction from A (q) C H"(q) to H"'(q) and 
satisfies len(y ? ) < x. 

Then, 

P u (A{j,o,m,x)) 

= P « U A ^b a > m,r o )oA 1 (j / 0 r m,r 1 ) o ••• o A q _i r Q _i)J. (5.11) 

' r 0' r l--- r Q-l e Q+ 


Since the paths y q are disjoint and piecewise linear, the conditions of the BK 
inequality are satisfied and we conclude that 

Q-i 

P u {A{j,a,m,x)) < e n P u(A q (j r a r m,r q )). 

r 0' r i' -' r Q-i e< Q+ ^=0 
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Noting that for a given configuration ( j,a / m ) the minimal length of a path 
connecting A(q) to H"(q) inside B(q) has the same distribution as Sm,N/ we 
conclude that 


Q-1 


Q-i 


P u (A(j,o,m,x)) < l n P u (X q (M,N) < r q ) = P,, Yj x ^ n ) < ■ 


r 0' r l'- ’ r Q-l e< Q+ 1?=0 
v’Q-i 

tf=0 r 1 <X 


,9=0 


which combined with | 5.10) , yields 

r. un t> \ \ „ V / number of choices 

P„(dist(0,Xei) < x) < > 


Q- M+N 1 


for (j, a, m) 


( Q-i _ 

<7=0 


< X 


(5.12) 


To complete the proof, we need to show that the number of ways to choose the 
triplets (j, a, m) is bounded by (2d(16M/N) d )®. To this end, assume that (j, a, m ) 
has already been chosen. By the definition of the times t"(q + 1), t q and 
f i;+ 1 and the points m(q) and m(q + 1 ), 


N\m(q + 1) - m(q)\oo < \Nm(q + 1) - y(t"(q + 1))U + | y(t"(q + 1)) - y(t q + i)U 
+ ly(Vi) _ y(tq)\oo + 1 y(t q ) - y(t"(q ))U 

+ | Nm(q) - y{t"(q)) U 

<3M + 5N. 


Since m(q + 1) e Z rf , it follows that there are at most ( 6^7 + lOj ways to choose 
m(q + 1) given m(q). Moreover, there are at most d choices for j(q + 1) and 2 
choice for a(q + 1 ), hence given ( j(q ), o(q), m(q)) there are at most 

/ M \ d / M\ d 

4n +1 °) s 4 6 n) 

choices for (j(q + 1 ), a(q + 1 ), m(q + 1 )) and in total at most 



choices for the whole sequence. 


5.3 Proof of Theorem 15.ll 

We separate the proof of Theorem |5.1| into two parts: a lower bound estimation 
and an upper bound estimation. We start with the first. 
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Formally, the claim for the lower bound is that for every u e [0, u t ) and every 
e > 0 there exists a constant c, depending only on d, u and e, such that for R 
large enough 

P„ (3* e ]R d : |*| = R, dist(0,*) < (q(n) - e)R) < e~ cR . (5.13) 


We start by showing that 11 ( 11 ) can also be obtained as the limiting distance 
between a point and a hyperplane: 


Lemma 5.7 For every u e [0 ,u,), 


dist(0 ,Hr) , s „ 

urn ---= tj(u), P u-a-s. 

R—> 00 J\ 


(5.14) 


Proof: From the definition of dist(0. Hr), 

dist(0, Hr) < dist(0,Rei), 


and therefore 

dist(0. Hr) dist(0,Kei) 

limsup--- < lim --- = i](u), P„-a.s. 

R->co K R ^°° K 

If 11 ( 11 ) = 0, then there is nothing left to prove since dist(0. Hr) > 0, hence assume 
that r\(u) > 0. Let co e Q be a realization such that 

lim inf ^— = i](u) - 25, for some 5 > 0. (5.15) 

R—>oo J\ 

Then, there exists an increasing sequence R; : that tends to infinity, a sequence 
of points zjc e Hr h and a sequence of paths y/ c from 0 to Zk such that 

dist(0,z,t) = len(yj;) < dist(0,HR t ) + 6 Rk < Rk(q(u) - 6). 

This however implies that 

dist(0,Zj;) R k (q(n)-5) 

limsup-—- < limsup-■— - < q(u) - 0 , 

k^>oo \ z k\ k^>0 O \ z k\ 

contradicting the uniform convergence of the distance function on S rf_1 proved 
in Theorem |4.4| Thus the event in 1 5.15) has probability zero and the claim 
follows. ■ 


Next, we show that the value of Sm,n defined in 15.41 cannot be much smaller 
than Mi](u). 


Lemma 5.8 For every e > 0 

lim rnaxP,, (s M n 5 M(q(u) - e)) = 0. (5.16) 

M->°o N<M 
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Proof: If i~i(u) = 0, then there is nothing to prove since sm,n > 0. We therefore 
assume that r\(u) > 0. Let 

i ' m.n = inf jlen(y) : y is a path from {0} X [0,N] d_1 to Hm} • (5-17) 

Since Tm,n is an infimum of path lengths over a set larger than that defining 
s m,n, it follows that tm,n ^ Sm,N/ therefore 

Pu(s m ,n < - £)) < P u (r M ,N < - e )) < P„(r M ,M < M(i](u) - e)), 

where in the last passage we used the fact that Vm,n is decreasing in N. Hence 
it is enough to prove that 

lim P„(r M ,M < - e)) = 0. 

M-»oo 


We now show that, in fact, it is sufficient to prove that for some fixed choice of 
5 = 5(e) e (0,1/2), 


lim P u (r M/25 M < M(ri(n) - e)) = 0. (5.18) 

M—>oo 


Indeed, let 5 e (0,1/2) and for v e Hq define vm(v) = dist (v,Hm)- Since {0} x 
[0,M] d_1 is contained in the union of boxes {{0} x (26 Mk + [0,26M] d_1 )); c eA(6) with 

A(6) = {k = (k 2 ,... ,k d ) e Z d “ 1 : 0 < fc, < ri/(26)l, V2 <i< d } 

it follows that 

rMM = inf(r M (i’) : v e {0} x [0,M] d_1 ) 

> min inf{r M (^) : v e |0) x (2 5Mk + [0,26M] d_1 )}. 

keA(5) 


Since the set A(5) is finite, and since each infimum for a fixed k has the same 
distribution as it is indeed sufficient to prove ( 5.18 1. 

Finally, let us prove < 5.18| . For every v e Ho such that < 2 5M we have 
dist(0, p) < 26M, and therefore 


r M (0) < r M {v) + 26M. 


Taking the infimum over all such p's yields 

Lvi(0) < £m,26m + 2 5M, 

and therefore for 5 < e/4, 

p u(r M , 2 bM < M(ri(n) - e)) < P„(r M (0) - 25M < Mffit) - e)) 

= P„(dist(0,H M ) - 26 M < M(t](u) - e)) 
<P u (dist(0,H M )<M(r ? ( M )-£/2)) 

which, by Lemma [5. 7\ tends to zero as M tends to infinity ■ 
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Proof of Theorem 5.1 (lower bound): For ;](n) = 0 the claim is trivial since 
P„(dist(0,x) < 0) = 0. For r](u) > 0, let N = min{M, [^JJ. For Q > - 1 

and f > 0, 


(Q -1 

Y J X q (M,N)<R(i 1 (u)-£) 

i=o 


= e /3R(r?(«)-£) . Ej( ^ 


-pXi(M,N) 


< e PR(i(“)-c) . Eu 
Q 


exp 


( Q-i 

-jS £x,(M,N) 

<fo 


< gP(Q+l)(M+N)(r;(H)-£) . ^ L r /JX 1 (M^)jQ 
= g^(M+N)(ij(«)-£) . ^(M+N)(rj(«)-£) . Ejj 

< eP(M + NMu)-e ). (gP(M+N)w»)-e). J e -fSM(7,(»)-e/2) + p a (x 1 (Ai,N) < M(j](u) - e/2)))) Q 

< JW+mm-e ). (e^ M£ /4 + g 2pMW»)-e)p II (x 1 (A^N) < M(ry(u) - e/2))) Q . 

(5.19) 


On the first line we used Markov's inequality; in the passage to the second 
line we used the fact that the XfM, N) are i.i.d.; in the passage to the third 
line we used the fact that R < (Q + 1)(M + N); the passage to the fourth line is 
an immediate algebraic identity; the passage to the fifth line follows from the 
inequality E u [e“^ x ] < e~^ a + P„(X < a), valid for every positive random variable 
X; finally, the passage to the sixth line follows from the choice of N. 

Recalling that Xi(M,N) ~ Sm,n, we obtain from Lemma |5.8| that 0 < P,,(X] (M, N) < 
M(i](u) - e/2)) < 1 for large enough values of M. For every such value of M one 
can find f3 m , depending only on M, d, e, u, such that 

(P„(X 1 (M / IV) < M(rj(:i)-e/2)))~ 1/3 < < (F^XfM^) < M(r](u)-e/2)))- 1/2 . 


When combined with 1 5.19) this implies 


(Q -1 

^ Xq(M,N) < R(i)(u) - e) 
q =0 


< £ MM+N)(tj(u)- f) 


• f(Pu(Xi(M,N) < M(r/(u) - e/2))) 24(rjfu)-e) + (P uiX^M'N) < M(jj(u) - e/ 2 ))) 1/2 j . 

(5.20) 


Using Lemma 5.8 one more time, we can choose M = M(e,d,u) large enough 
so that 

(P„(Xi(M,N) < M(i](u) - e/2))) W + (P u (Xi(M,N) < M(tj(w) - e/2))) 1/2 


< 32 d ■ max < 2, 


&?(») j 


r-Kr- 


(5.21) 
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For such choices of M and j3, we get from 1 5.20\ , 1 5.21} and Proposition|5.6| 


d\Q 


P„ (dist(0, Re{) < - e)) < e P( M+N M u )-e) . L M^rKr 

o>—^_1 ' ' 

< e m+mnw-e ). ^ 2~ dQ 

Q^h+n - -^ 

< e m+mnw-e) . 2-i+d+i, 

Recalling that M, N and [’> are fixed, this gives the desired exponential decay in 
R. 

Given the result for x = Re i we turn to deal with general points, x e R ri , \x\ = R. 
Due to the invariance of P u under rotations, dist(0, x) for |x| = R has the same 
distribution as dist(0. Re i), and therefore for large enough R 

P u (dist(0,x) < (/](«) - e/2)R) < e~ c ' R , Vx e JR d , such that |x| = R. 


Taking an e/2-net At on S d 1 such that |Af| < we get that 

P„(3x e N such that dist(0,Rx) < (n(«) - e/2 )R) < e~ ClR . 

e d 

For every x e ]R rf such that |x| = R there exists a y e At, such that \x/R-y\ < e/2, 
and therefore 

dist(0, Ry) < dist(0,x) + |Ry — x\< dist(0,x) + ^R. 

Hence, 

P„ (dx e lR rf : |x| = R, dist(0,x) < (/](») - e)R) 

< P„ (3 y e At such that dist(0, Ry) < ( ij(n ) - e/2)R) 

< - d e- c ' R , 

£ d 

which concludes the proof. ■ 


Next, we turn to prove the upper bound in Theorem 5.1 which states that for 
every u e [0, u t ) and every e > 0 there exists a constant ci, depending only on 
d, u and e, such that for R large enough 


P u (3x e R‘ ( : |x| = R, dist(0, x) > (rj(u) + e)Rj 


< e 


,-ciR 


(5.22) 


We start with the following lemma. 
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Lemma 5.9 For S < Rlet 


. , f, , , y is a path from Se\ to Rei 
Ps, R = mf jlen(y) : such that Hg K y K Hr 


Then 


r E, ( [ Po ,r] , . 

p lim —n— = '?(“)■ 

R—>oo K 


Proof: We follow [ HW65 . Theorem 4.3.7]. Since dist(0, Rt’i) < po,s 
liminfhfeM! > lim E„[dist<0.K ei )] = 

R— »oo R R—> DO R 


For A > 0 define 


,, . , f, , , r is a path from Sei to Rei 

Ps,R = mf -j len(y) : such £ at Hs _ fc < y < Hj 




R+i 


(5.23) 


so in particular, p$,R = Psr- Since for every A > 0 and 0 < R] < Rt, 

Po,Ri+R 2 — Po,Ri + PRi,Ri+R 2 

one can apply the subadditive ergodic theorem lKin73 l for p k 0 R to obtain 


(5.24) 


E M,r\ 


lim 

R—>oo K 


= rf{u) 


(5.25) 


which, by the definition of p k Q R/ satisfies 

Tj(n) < rf(u) < rf~ l (u) < T]°(u), VA > 1. 

Noting that for every R, k > 0 

P-k,R+k 5 p-kfl + Po,R + PR,R+k/ 

it follows by fixing A, taking expectation, dividing by R and taking the limit 
R —> oo that 

'7°(H) < 

which together with | 5.25] implies that t] 0 («) = if (it) for every fixed A > 0. Since 
for every fixed ae Q and R > 0, p 1 ^ R is a monotonically decreasing function in 
A converging to dist(0, Rt’i), it follows from the monotone convergence theorem 
that 

lim E„[po R ] = E„[dist(0,Rei)]. (5.26) 

k —>oo ' 

Using once again the sub-additivity ||5.24] we conclude that 


E uiP\ 


0,R J 


R 


> rf(u) = if(u) 
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for every R > 0 by Fekete's subadditive lemma. Combined with ( 5.26) this 
yields 


E»[dist(0 / .Rei)] 0 

R — « v /* 


Letting R 
proof. 


implies rf(u) < r\(u). Together with 15.231, this completes the 


Proof of Theorem 5.1 (upper bound): Proving the upper bound is in fact much 
simpler than proving the lower bound. Fix e > 0. By Lemma 5.9 there exists a 
sufficiently large Rq > 2 such that 


E [Po,J , , £ 


(5.27) 


For i > 0 let X, = p/R 0 ,(i+i)R 0 - The random variables (X;),>o are i.i.d. (with the 
same distribution as po,R 0 )/ and 

LR/RoJ-l 

dist(0,Rei)< ^ X,+R 0 . (5.28) 


i =0 


Using (]5.2 71 and (5.28) we deduce that for every R > 5Ro/s, 


fLR/RoJ-l 


<P U 


P„(dist(0 / Xei) > (rj(w) + e)R) < pj ^ X, > (rj(u) + e - R 0 /R)R 

i =0 ) 

b RrjJ-1 „ c r v 

£ (Xi - E„[x,]) > (,](«) + £- Y - -4^) 

bRoJ-l ~ 

£ (X, - E„[X,]) > - ( <;R 


R 


<P„ 


V i =0 


Applying the function x i—> eJ ix (for some jS > 0) to both sides in the last term, 
using the Markov inequality and then the independence of the X,'s, we can 
bound the last term on the right-hand side by 




e 5 






= e 5 


E u [‘ 


,^(PO,R 0 -E„[PO,R 0 ]) L R ° 


^(POJ! 0 - E »[PO,R, 


1 

k,D]*o. 


(5.29) 


<e~f R -E u [e 

Since E„ ^e/ J lf'o,Ro _E »[po,R 0 l)j < £ J( ^P°. R 0 J < eP R ° < oo we can expand the exponential 
inside the expectation into a power series in f and obtain that 

E„ [ e / s (P°.R 0 - E dpo,R 0 ])j _ I + 0(13 2 ). 

Hence ( |5.29) is bounded by exp - -y- jxj- By taking j 3 = f(e) > 0 small 

enough we can make the last term exponentially decaying in R, thus completing 
the proof. ■ 
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5.4 Proof of Theorem 15.21 


The proof of Theorem 5.2 also follows ideas of Kesten ! Kes86 1. Unlike the proof 
of Theorem 5.1, some additional elements are needed in order to apply to the 
continuous case. We start with some preliminary results. 

Lemma 5.10 There exists a constant C 4 = Ci(d) such that for every M > 1, the 
Euclidean length of every true geodesic contained in a box of side length M is at most 
CiMf 


Proof: Denote the box by Bm ■ Let a> e LI be a configuration, and let y be a true 
geodesic contained in Bm■ Since lcn(y) < diam c (Bjvr) = VrfM, the Euclidean 
length that y acquires in Bm \ <S(a>) is at most VrfM. Therefore, it suffices to 
show that the Euclidean length that y acquires in S(co) Pi Bm is 0(M rf ). 

Let {A,-}^ be the connected components of S(a>) P Bm which y intersects. Note 
that y P A; is connected (otherwise it would not be a true geodesic), so there are 
well-defined entry and exit points to y P A;. Denote by xl e A, (resp. x^ e A,) 
the point in supp(<x>), from which the entry point (resp. exit point) of y P A, is 
of distance at most 1. Let xl = xj,xf,..., x k . ! = xb be points in supp(ai) P A„ such 
that \xl - x^\ <2 if \k - j\ < 1, \xl - x^\ > 2 if \k - j\ > 1. We can always find such 

points by taking supp(cc) P A, and omitting points. By construction, xl e A, for 
every j, and therefore the length y acquires in A, is at most 2k, (the length of 
connecting xl with xl +1 plus the length of connecting the entry and exit points 
with x 1 and x L ). 

The length y acquires in S(co) P Bm is therefore 2k\ + ... + 2/c ; \;. However, by 
construction, the unit balls centered at (xy , );</v / 2 ;<) Cl are mutually disjoint. By 
a volume consideration, the number of disjoint unit balls in Bm is at most 
K d M d , where K d is the volume of the Euclidean unit ball in !R". Therefore 
k\ +... + kjv < 2 K d M d , and so the total length y acquires in S(co) P Bm is bounded 
by4KrfM rf . ■ 

We also recall the following variant of a result by Roy |Roy90J ; see also 1MR961 
Lemma 3.3]: 


Theorem 5.11 ('foydO/) For x e ]R‘ f denote 


S(x;co) = 


the connected component ofx in S(co), 

0 , 


x e S(co) 
x t S(co) 


Then, there exist for every u < iu positive constants C,c depending only on u and d 
such that 


P„(diam c (vS(0;cL>)) > t) < Ce ct . 


In order to prove Theorem 5.2 as well as for future use, we will need a stronger 
version of Theorem l5.11l 
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Lemma 5.12 Let 

W= |J 5(x ;( y). 

xe[0 ,lf 

Then, for every u < n, there exist positive constants C$,C(, depending only on u and d, 
such that 

P„(diam c (W) > t) < C^e~ Cbt . 


Proof: Observe that every unit ball B(x, 1) for x e R d contains a point in the grid 
. Therefore, each connected component in W must contain a point from the 
set 

L = [q 6 : \q\ 2 < 2 } 

whose size is finite (and depends only on d). 

Assume that diam c (W) > t. Then, there exists a connected component whose 
Euclidean diameter is at least (f /2 - Vrf). By Theorem 

P u (diam e (W) > t) < ^ P u (diam e (<S(x)) > f /2 - Vd) < Ce~ d . 

xeL 


5.11 


Proof of Theorem [T2j The proof follows by an argument very similar to the one 
used in the proof of Theorem 5.1 We start with Part (1), showing that there 
exist constants a, C 2 and c 3 , such that for every R > 0, 


3 a true geodesic path starting at 0 
such that len e (y) > R and len(y) < ^ len c (y) 


< C 2 e 


-c 3 R 


First observe that 

P 


3 a true geodesic path starting at 0 
such that len c (}') > R and len(y) < - len c (y) 


S=LRJ 


3 a true geodesic path starting at 0 
such that len e (y) e [S, S + 1] and len(y) < yf S + 1) 


so it is enough to show that there exist constants a, C, C 3 such that for every 
natural number S 


3 a true geodesic path starting at 0 
such that len c (y) e [S, S + 1] and len(y) < -kS + 1) 


< Ce 


-C3S 


Fix S e N and let y be a true geodesic path starting at the origin such that 
len e (y) e [S, S + 1]. Fix M e IN and N = 1 and define the sequences f, ; and x I? , 
q = 1 ,..., Q in the same way as in the proof of Proposition|5.6| By construction. 
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each of the segments y|[t,,t ?+1 ] is contained in a box of side length M + 1. By 
Lemma 5.10 the Euclidean length of a true geodesic in each of the boxes is at 
most C 4 (M + l) d , hence 

Q 

Q > -- - 1. 

C 4 (M + l) d 

By repeating the argument of Proposition |5.6| we get that for every f} > 0 

3 a true geodesic path starting from 0 
such that len e (y) e [S, S + 1] and len(y) < ^(S + 1) 

o f Q_1 1 

< Y, (2d(16M) d y~ P u Y X ^ M ’ 1) < -( S + 1) 


Q> 


C 4 (M+1)“ 


7-1 


,7=0 


< Y (2d(16M) d )V (S+1)/a E„ 

\Q 


Q-i 

n- 

q=0 




= g/3(S+l)/a 


^ (2d (16M) d ) E„ [^ X,W1) ] . 


Q> 


C 4 (M+1)“ 


7-1 


With W as in Lemma [5.12[ 

P„ |diam c (W) > yj < C 5 e 


-c 6 M/2 


Each connected component of S is composed P„-a.s. of finitely many balls; see 
Lemma 4.1 Hence connected components of S, as well as W, are compact sets. 
Therefore, one can find a small enough 6 = <5(M) > 0 such that 

P„(dist(W„S \ W) > 6(M)) > 1 - C 5 e~ CbM/2 . 

It follows that with P„-probability at least 1 - 2C^e^ Cf ' M ' 2 the set W is contained 
within the box [-M/2,M/2+l] d and is at a distance at least 5(M) from any other 
cluster of S, which in particular implies that e pXi(M,i) < e fib. Consequently, 

+ 2C 5 e~ C(,MI2 


and thus 


3 a true geodesic path starting at 0 
such that len c (y) e [S, S + 1] and len(y) < ^(S + 1) 

< gPiS+iya 

Q> ; . i 


(2d(16M) d ) Q (e-WW + 2C 5 e- C6Mlz f 


“ C 4 (M+ l) d 


TakingM large enough so that 2d(16M) d -4Cse '■ v,iV1/2 < |, and then /I large enough 
so that e~^ M ) < 2C^e~ Cf,MI2 , and finally a large enough so that e^ a < 2 c * (M+1)d 
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gives 


3 a true geodesic path starting from 0 
such that len e (y) e [S, S + 1] and len(y) < -(S + 1) 


J < (1) C4( " +1) " 


proving Part (1) of Theorem 5.2 


We next prove Part (2), namely that for every x, y £ R d and for every R > a\x — y\ 
P„ (3 y e T 0 (x,y) : len c (y) > R) < C 2 e~ C3R . 


By the invariance of the measure P„ under translations and rotations, it is 
sufficient to prove the result for the case x = 0, y = Se 2 for some S > 0. If 
y £ To(0, Se i) is a path such that len 0 (v) > R > aS, then we can find times f, ; and 
points x q ,q = 0,...,Q, as in Proposition |5 .6 1 (vvith N = l),withQ > aS(M + l)-l. 
Since len(y) = dist(0, Se i) < S the result follows by the same argument used to 
prove Part (1). ■ 


6 Further result for uniform point distributions 


In this section we exploit the results obtained in the previous two sections to 
prove more results for uniform point distributions. These include the geometric 
concentration of geodesics and further properties of the function 


6.1 Geometry of geodesics 

In order to obtain results on the limiting distance in the setting of general 
intensity measures, we will need control over the geometry of geodesics, or 
approximate geodesics, connecting pairs of points. For x, y £ 1R 1? , we call a path 
y between x and y an e-geodesic if 

len(y) - dist(x, y) < e\x - y\. (6.1) 

We denote by T c (x, y) the set of e-geodesics which are also Euclidean geodesics 
inside S. In particular. To (x,y) is the set of true geodesics between x and y 
defined in Section[5] 

The goal of this subsection is to show that there exist geodesics that do not 
deviate significantly from Euclidean segments: 

iProposition 6.1 Let n £ [0, «*). There exist C 7 ,cs > 0 such that for every e > 0 and 
sufficiently large (depending on e) \x - y\, 

P„ (Vy e T £ (x, y) : d H (y,[x,y]) > e\x - y\) < C 7 e~ c * ]x ~ y]V2 , (6.2) 

where [x, y\ is the linear segment connecting x and y. 
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Proof: We follow the proof of IBI .PR 1 5 Proposition 3.2]. Fix e > 0 and x, y e R rf . 
We will show that if \x - y| is large enough, then there exists with P„-probability 
> 1 — a curve y e T f (x, y) satisfying 

dn(y,[x,y]) < e\x-y\. 


Let N = F20a/el, with a as in Theorem |5.2| Define the sequence of vertices z k , 
k = 0,..., N, along [x, y] by 

Zk = ( 1 ~^) X+ ^ y ' V0 ^^ M 

For 0 < * < N - 1, let y* = y k (co) e To(z k ,z k+ \) be a true geodesic connecting 
Z/c and Zk+i and define y = (■■■ (yo * yi) *•••*) * }'n-i) to be their concatenation, 
connecting x and y. 


Since \z k - z k+k \ = \x - y\/N is of order \x - y\ it follows from Theorem 5.1 that 
once \x - y\ is large enough 

\dist(z k ,z k+1 )-Tj{n)-\z k -z k+1 \\<^\z k -z k+1 \, VO < fc < N - 1, (6.3) 

with P„-probability at least 1 - Ne -Cll * -yl . Similarly, whenever \x - y\ is large 
enough we have with P„-probability at least 1 - that 


| dist(x, y) - i](u) ■ \x - y\\ < -|x - y\. 
Consequently, under the events in (|6.3|l and ]6.4), 


(6.4) 


N-l 


N-l 


len(y) < E len(yj.) = dist(z fc ,z fc+1 ) 

k=0 k=0 

N-l 

< ^ yi(u) • \Zk - z k+ i\ + -\z k - Zfc+ilJ 


(6.5) 


k=0 


= '?(») ■ \x - y I + ^\x- y\ 
< dist(x, y) + e\x - y\, 


which implies that y e T. (x, y). Finally, by Theorem 5.2 2) we have with P„- 
probability at least assuming |x-y| is sufficiently 

large, that 

( 6 . 6 ) 


leneOht) < — |x-y|, W0<k<N-l 


which implies that for every 0 < k < N — 1, 

y k c [Zjt,z i+ i] + B ^0, ~ \x - y|j. 


and thus 


y c [x,y]+B^0, ||x-y|), 
thus proving the existence of the required e-geodesic. 


(6.7) 

( 6 . 8 ) 
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6.2 An upper bound on r\(u) 

The goal of this subsection is to prove an upper bound on the function !](n): 
‘Proposition 6.2 For every u e [0, u») 

i](u) < e~ UKi . 


Proof: First observe that x e S(o>) if and only if there is a point y e supp(a)) n 
B(x, 1). Since the the number of points in the support of co inside B(x, 1) is 
distributed like a Poisson random variable with parameter UKd it follows that 


P u (x g S{co)) = e~ UKi . 


For R > 0 let y R : [0,1] —> IR d denote the path yi<(t) = tRe i. Then 


dist(0,Re!) < len(y R ) 


i r R 

~ p J Leb [ (dx). 


R ~ R 

Taking expectation on both sides and using Fubini's theorem, 
E„[dist(0,Rei)] 


R 


i r 

e -UK 

R Jo 


UK * Lebi (dx) = e~ UK “. 


Taking the limit R —» oo and using Theorem [44] gives r](u) < e UKd . 


tftemarlf In fact, one can show that the bound e~ UKd is not tight, and that 
/](«) < e~ UKd actually holds, by considering several paths instead of one as 
follows. Take large balls around 0 and Re\ (whose radii are independent of R). 
An < -bound is then obtained by considering several paths from 0 to Re 
which are at distance > 2 from each other outside these balls. 


6.3 Continuity of i](u) 

Proposition 6.3 The function tj : [0,«,) —* (0,1] is continuous. In addition, i] is 
monotonically decreasing in u and ?](0) = 1. 

We start the proof by introducing a natural coupling of the probability measures 
P„ for u > 0. Let 

~ _ y Xi e ]R rf , U{ e [0, oo) for all i > 0 and a>(A x [0, u\) < oo ) 

I 2-j ‘ for all bounded, Borel-measurable A c ]R‘ < and u > 0 J 
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and for u > 0 and co e Q define 

<0 U = 5 Xi . 

(Xi, Uj ) € Slippy 
Uj < U 

Note that co u e Q, so dist(-; co u ) is well-defined. Finally, let P be the probability 
measure on Q under which a> is distributed like a Poisson point process on 
]R rf x [0, oo) with intensity Leb n(dx) X Lebi(dx)|[o,oo). 

One can verify that for every u > 0, the distribution of co under P„ is the same 
as the distribution of co u under P. That is, we constructed a coupling of the 
probability measures P„ for u > 0 under which 

a) u < co U ', VO < u < u'. (6.9) 

We will denote by E expectation with respect to the probability measure P. 


Proof of Proposition 6.3 Fix u e 
[u - 5o, u + <5o] c [ 0 , iu) if u > 0 

compact interval by I. 


[0, !(*) and £ > 0. Choose <5o > 0 such that 
and [0,5 q] c [0, iu) if n = 0. We denote this 


Denote 


A r = {3y 6 r 0 (0,Rei) : len c (y) > aR} U j| dist(0,Pei) - i](u)R\ > eR/6}, (6.10) 
where a is as in Theorem |5.2[ 2). 

For every u' e f, using the fact that i](u), i)(n') e [0,1], we have 


\r](ii') ~ ri(n)\ < |e [(jj(m') - r](u))l^(w. i ')lA^(««)]| + P (A R (co u >)) + P(A r (oj„)) 
E [[dist( 0 , Rei; a> U ') - dist( 0 , Rep Si,)] • 1 a^(«»') 1 a^(S'„)]| 


1 

< — 
“ R 


+ 3 + + Y(Ar(co u )). 


( 6 . 11 ) 


By going back to the proofs of Theorem |5. 1 1 and Theorem |5.2[ 2) one can verify 
that both a and the constants C\, C2, C3 can be chosen uniformly on the compact 
interval I. That is, there exist positive constants a, Cg and C 10 depending on u 
and 60 such that for every u' e I, P„'(Ar) < Cge~ CwR . 


Due to the uniform bound on the probability of Ar on the interval l one can 
choose R large enough (depending only on £ > 0 and I) so that 

P(A*(2v)) + POMS*)) < ( 6 . 12 ) 

I, it is enough to show the existence of 5 > 0 such 

that for \n — u'\ < 6 , 


Combining (| 6 . 11) and 16.12 


1 

R 


E [[dist(0,i?ei; av) - dist(0, Rev,co u )] ■ 1a'(w«')1a|(w u )] 
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This holds since on the event a) U ’,oo u e A C R the values of dist(0, Re^;a> u f and 
dist(0, Rc’i; co u ) are the same whenever 

supp(Sv) n B(0,Ra + 1) = supp(oi u ) n B(0 ,Ra + 1), 


and therefore 


1 

R 


E [[dist(0, Rei; ov) - dist(0,Rei;cu H )] ■ 1 a c r (co U ’)1a c r (co u )\ 


< 2 ■ P(supp(av) n B(0,Ra + 1) + supp(ay ( ) D B(0,Ra + 1)). 


Assume without loss of generality that u' < u. Due to the coupling, the point 
measure co u is obtained from 2v by adding to it an additional independent 
point measure in Q which is distributed as a Poisson point process with intensity 
measure (u-u')-Lebrf(dx). In particular, the probability that there is an additional 
point inside the ball B(0 ,Ra + 1) is 


l _ e -(u-u')K d (Ra+If < | _ e -6K d (Ra+l) 1 ' 

Thus, with probability at most 1 - e~ 6Kd( - Ra+1)J the point measures to U ’ and a> u do 
not coincide inside the ball B (0, Ra +1). Recalling that R depends only on e and 
I we can choose 5 small enough so that 2(1 - e~ 6Kd ^ Ra+ A li ) < e/3, thus completing 
the proof. 


6.4 Volume convergence 

Proposition 6.4 For every u > 0, P u -almost surely 


lim v K ([0,1]“) = lim 

R —>oo M—>oc 


Lebrf([0,M] d \ S(a >)) 
Lebrf([0,M] d ) 


= e 


= ^a(H)([0,l] d )- 


(6.13) 


Proof: As observed in Proposition |6.2| 

P„(0 £ S{co)) = e~ UKd . 


Using the ergodicity of the model, see Lemma 4.5 and the ergodic theorem we 
can conclude that 


lim UMI0,Mr\.SM) = 

m-> o° Lebd([0,M] d ) oo Lebrf(B(0, 

= E t ,[l 0 ^)] = e~ UKd ■ 


3(0 ,M))X M 1 " 5< " |Leb ' ( * ) 
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7 Convergence for uniform point distributions 


In this section, we exploit the results proved in Sections 4j|6]to prove Parts 1 
and 4 of Theorem 12.21 for the case when u is uniform and D is convex. These 
assumptions are relaxed in the next section. Since for uniform u, we have a 
natural coupling of the measures P u ,r using a single measure P„, we will have 


in fact a slightly stronger result than stated in Theorem 2.2 
The main result of this section is the following: 


‘Theorem 7.1 Let u e [0, iu) and let D c R d be a convex, compact d-dimensional 
manifold zvith corners. Then 


x,yeD 


lim sup d°(x,y) - dist%(n R (x),n R (y)) 

<— »CX> _7-V ,V ' 


= 0 , 


P u-a.s. 


(7.1) 


In particidar, since n R is onto, the sequence (D R , dist°) converges P u -a.s. to (D, 
with respect to the Gromov-Hausdorff metric, where shoidd be interpreted as in 
|2T3). 


Proof: Denote = (D R , distg), = (D K ,dist R ) and C D = (D,d° () ). We 
denote by the projection D —» D R when considered as a mapping between 
C D and Jdj). We will use the space B^ as an intermediate metric space in order 
to bound the distortion of n l f We denote by fi 1 ^ the projection D — > D R when 
considered as a mapping between C D and B^l and by Id the identity from D R 
to itself when considered as a map between B 1 ^ and TPff. See Figure |3] for an 
illustration. 


C D = (D,d° u) ) ---= (D*,dist°) 



= (D r , dist R ) 

Figure 3: The spaces and C D . 


By the triangle inequality, 

dis 7i^ < dis + dis Id, 


where the distortion is defined as in Subsection 2.3.6 Thus, it is enough to 


prove that with P„ -probability one, both dis fi^ and dis Id go to zero as R —» oo. 


For dis n^, observe that D is compact and convex, and since u is constant, it 
follows that D is convex with respect to d n ( U ), hence The mapping 
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4.2 


has P„-a.s. an asymptotically 


: Cf'j —* is onto and by Theorem 
vanishing distortion. 

The rest of the proof shows that the distortion of Id : - » also vanishes 

asymptotically. We use the concentration result Proposition [64] as follows: Fix 
5 > 0, and denote by D(6) the set of points in D whose distance from the 
boundary 3D is greater than 6. For R > 0, let Nr^ denote a finite 1 / VR-net 
of D(5) such that \Nr^\ < CR d/ 2 and for every x, y e Nr#, \x — y\ > c/ VR, with 
C and c depending only on d and D. Denote 0 < £ = 6/(2 diam(D)), and let 
x, y e N Rj6 . By Proposi Mon |6. 1 1 we have that for R > Ro(6), 

P» (Vy e T c (Rx, Ry) : d H (y, [Rx,Ry]) > eR|x - y\) < C 7 <r C8(R| *^ l)1/2 

< Cge~ CwRlli , 

and therefore by a union bound argument and the fact that \Nr ^l 2 < C R lf , we 
have 

P„ (3x ,y e Nr s Vy e T £ (Rx,Ry) : d H (y, [Rx,Ry]) > eR|x - y\) < CnR rf e _Cl ° Rl/4 . 

(7.3) 

Considering the sequence of events in ( ]7.3) with R replaced by m e N, we get 
that the sum of the probabilities is finite, hence by the Borel-Cantelli lemma, 
we have that P„-a.s. there exists Mo(6, a>) such that for every m > Mo(6, co), 

mb 

Vx, y e N m £ By e T £ (mx, nzy) : d H (y, [mx, my]) < em|x - y| < —. 

Since D is convex, [x, y] c D. Therefore by the definition of 6, it follows that for 
every m > Mo and every x, y e N m ,g there exists an e-geodesic with respect to 
dist R that remains in Dr for R large enough, hence 

lim sup | dist m (x,y) - dist°(x,y)| < e, P„-a.s. 

m—>oo _■» T 

Since N m ^ is a 1/ yfm -net of D(6) it follows that 

lim sup | dist m (x,y) - dist°(x,y)| < e, P„-a.s. 

m —.too 

x,yeD(6) 


2 diam(D ))' 


Ptra.s. 


and therefore 

lim sup | dist m (x, y) - dist°(x, y)|<e+45 = 64 + 

m—>DO x,yeD \ 

Since for every R > 0 there exists an m such that \m — R\ < 1, it follows that for 
such a choice of m, for every x, y e D, 

dist(mx, my) dist(Rx, Ry) 


| dist m (x, y) - dist R (x, y)| = 


m 


dist(mx,Rx) dist(my,Ry) 

<-h-h dist(Rx, Ry) 


m 


|x| D 


lylo 


m 

h - ylo 


R 

1 _ 1 
m R 


R - 1 R - 1 R - 1 


4su P ze p I 2 1 

R - 1 
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Similarly, 


|dist“(x, y) - dist“(x,y)| < S ^ ■ 

It follows that 

limsup sup |dist R (x, y) - dist°(x, y)| < e + 46 = 6(4 + ^ |, P„-a.s. 

R->co x,yeD \ 2diam(D)/ 

Since 6 is arbitrary we finally have that 

lim sup | disbi(x, y) - distj^x, y)| = 0, P u -a.s. (7.4) 

R ^°°x,yeD 

This shows that the identity mapping Id : S R —» Jl R has an asymptotically 
vanishing distortion P„-a.s., which completes the proof. ■ 

The following proposition proves Part 4 of Theorem |2.2| for the case of constant 
u, that is, that (n^) -1 : D' R —» D is asymptotically surjective. 

Proposition 7.2 

lim d H (D, n K \D' K j) = 0, P u -«.s. 

R—>oo 


Proof: It follows from Lemma 5.12 that 


I 3a connected component C in U.re[-M,My d>(x; co) 
“ ( such that diam e (C) > log 2 M 


* E '■ 

zeZ li n[-M-l,M] li 


3 a connected component C in (Jjez+miy *S(x; co) 
such that diam c (C) > log 2 M 


<C 5 (2M + 2) d e- C6lo s 2(M) . 


The choice of log 2 M is dictated by the need of satisfying two conditions: we 
need a term which is o(M), and the probability must decay sufficiently fast; see 

0 

Taking M = KR, with K = diam t .(D) and R > 0, we get that with P„-probability 
at least 1 - Cf 2KR + l) rf e —C6lo g the distance from any point in S IT R ■ D to 
R ■ D \ S is at most log 2 (KR). Thus 


P u \d H (D,(n^-fD' R ))> 


log 2 (^) 

R 


< C 5 (2 KR + 2) rf e - c <dog 2 (iad 


(7.5) 


Given e > 0, consider the sequence (dn(D , (n° f ) 1 (DJ n£ ))^ . Since the right- 

hand side of (73} is summable in m, it follows from the Borel-Cantelli lemma. 
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for all but finitely many 


that P„-almost surely d H (D, (n° E ) 1 (D' m£ )) < 1 ° s £ ( f £ ”' ) 
m's and in particular that 

lim d H (D, (tt° £ ) -1 (D; )) = 0. 

m—>oo 


Since for every 0 < R < S we have d H ((n R )~ 1 (D' R ),(n 1 ^)~ 1 (D' s )) < [1? - S| and 
since the sequence (dn(D, (n^)' (D' mi .)))^ is e-dense in [0, oo) it follows that 
P„-a.s. limsupp^^, c1h(D, (7 i?) -1 (D')) < e. Since e > 0 is arbitrary we get that 
P u -a.s. lim R ^ d H (D, (n°)-\D' R )) = 0. ■ 


8 Proof of Theorem 12.21 


In this secti on w e prove Parts 1,2 and 4 of Theorem 2.2 Part 3 was proved in 
Proposition |6.2| Let D C R rf a compact d-dimensional manifold with corners, 
and let it : D —> [0,«.) be a continuous function. Since D is compact, we 
can always extend it continuously to IR'^ without enlarging its upper bound. 
Therefore, we can assume without loss of generality that u : ] R rf —> [0,»,), with 
sup u < u t . The parameters u and R > 0 define a process with probability 
measure P„ r. Similarly to the proof of Proposition |6.3| we start by introducing 
a natural coupling of the probability measures P„,r for a given R > 0. 


Let 


Id — i CO — &{xi,ui) ■ 


i> 0 


Xi e R rf , Ui e [0, oo) for all i > 0 and co(A x [0, u]) < 
for all compact A c R rf and u > 0 


and for a continuous function u : R ri 

C0 U = 


-* [0, oo) and a> e Q define 

E 


(xj, Uj) e supptd, 

Ui < ll(Xj) 


Finally, let P R be the probability measure on Q under which a) is distributed like 
a Poisson point process on R rf x [0, oo) with intensity R d Leb a{dx) x Lebi(dx)|[o,a>). 

One can now verify, that for every continuous u : 1R‘ ? —> [0, oo), the distribution 
of co under P UjR is the same as the distribution of a) u under P R . That is, we 
constructed a coupling of the probability measures P u R , for u : R rf —» [0, oo), 
under which 

co u < a>w, Vh,h' : R rf —> [0, oo), u < u'. 


8.1 Metric convergence 


In this section we prove Part 1 of Theorem 2.2 which with the coupling con¬ 
structed above states that: 


50 







For every e > 0, 


lim Pr 

R —>oo 


sup 

X rV €D 


\tfo u (x,y) ■ 


dist°(7T R (x),n K (i/);a;„) < e 


= 1 . ( 8 . 1 ) 


In particular, this implies that for every e > 0, 

lim P R (dGH({DR,distz(-;aj u )),(D r d°)) > e) = 0. 

R—>oo 1 

Theorem |7.1| states that this holds when D is convex and u = i/q is uniform. 

We start with a proposition showing that if D is convex and u : D — » [0, u t ) 
only takes values within a small interval [u m in, »max]/ then we can bound 
(with high probability) the distortion of the projection n R between (D,d® ou ) 
and (D R (co u ), dist R (-; at,,)). 

“Proposition 8.1 Assume that D is a compact, path-connected and convex set with a 
non-empty interior, and let u m \ n = mino u > 0 and u max = maxo u < u,. Consider 
n R as a function (D,d° ou ) -> (D K (S' u ),dist“(-;aT„)). Then 

lim P R (dis tir > 3Ai] • diam c D) = 0, 

R—>oo 

where At] = t](u m „) - i]{u max ). 


Proof : Consider the coupling as above between co u , co Umin and co Umax , with u m i n 
and « max viewed as constant functions. Then, 


di st R (-;<W m J > dist°(-; w u ) > dist£(-; co Uim f) and d ° (Umin) > d l q { u > d“ (lw) , 

and in particular, 

| dist R (-; o; Um J - dist R (-; a;„)| < | dist°(-; w Uw J - dist R (-; w Um J\, (8.2) 

and 


I - <o J < I< (1W) - < ( , w) 


< At] ■ diam c D. 


(8.3) 
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Using l |8.2[ and ( |8.3j i, and by a repeated application of the triangle inequality, 
dis n R = sup | d° ou (x,y) - distR(x,y;a?„)| 

x, v eD 

< sup | d° ou (x,y) - d° Umjn) (x,y)\ + I- dist^(x,y;ST Umin )| 

x,yeD 

+ | dist°(x, y; af Umln ) - dist°(x, y; co u )\ 

^ SU P - d °(u m J x ' l j)l + I d n(u m J x 'V) ~ dist°(x,y;a; Umin )| 

x,yeD 

+ \ distjfo, y; af Umln ) - dist£(x, y; S„ max )| 

^ SU P 21 d° (Umax) (x, y) - d° (Umin) (x, y)| + 2|d° (Umiii) (x, y) - distjfa, y; S?„ mta )| 

x,yeD 

+ I d n("™A x ' V) - dist°(x, y; S,, max )| 

< 2Aq ■ diam,. D + 2|d° Umjn) (x, y) - dist°(x, y; Sf Umin )| 

+ \ d n(u m „){x, V) - dist°(x, y; S' 1(max )|. 

(8.4) 


The result now follows by applying Theorem 7.1 to the last two addends. 


We now prove Part 1 of Theorem 2.2 The idea behind the proof is the fol¬ 


lowing: We partition D into small convex sets, in each u varies only a little. 
Proposition |8. 1 | states that with high probability, the distortion of the projection 
in each set is small. We then glue the sets together and show that the accumu¬ 
lated distortion remains small. A technical complication arises when D cannot 
be partitioned into finitely many convex sets. We overcome this problem by 
considering sets slightly larger and slightly smaller than D, denoted by D and 
D, that can be partitioned in such a way. 


Step I: Partitioning D. Let n be a large natural number to be chosen later 
(independent of R). Cover D with cubes of edge length 1/n, whose corners 
are on the lattice henceforth, "vertices" refers to the corners of the cubes. 
Denote the cubes that intersect D but not <9D by a n ,\, • • •, □«,£„, an d those that 
intersect dD by □„ *.„+!,..., a„^ n+mn . Since D is compact, there exists a slightly 
larger compact set □ containing 1J ."| m " for all n > 1. The function /] o u 
varies on each cube by some A n j. Denote A„ = max, A ;1/ ,. Since tjo u is 
continuous (Proposition |6 . 3) and □ is compact, A„ —> 0 as n —> oo. Denote by 
T„ the union of the facets of the cubes and let 

k n k n +m n 

D = |^J ^n,i/ D — ^n,i‘ 

i =1 i= 1 
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It follows from the definition of D and D that DcDcD and both | diam D - 
diamD| and | diam D - diamD| are of order 1/n. Moreover, we claim that 

iim sup \d° ou (x, y) - d° ou (x, y)\ = 0. (8.5) 

n ~‘*°° x,yeD 

Indeed, note that d[! u < dy au for any pair of points in D, and therefore we only 

need to prove that for every x, y e D and every simple curve y c D between x 
and y, there exists a simple curve y' c D between x and y such that 

len d,JV) < len dqou (y) + o(l), (8.6) 

where o(l) is with respect to n, independent of R, x and y. To prove ( ]8.6) , 
note that since D C lR rf is compact, D and D are homotopic whenever n is 
large enough ( dD is a compact submanifold with comers and therefore its 1/n- 
neighborhood is homotopic to itself for n large enough) and d//(D, D) < xfd/n. 
Then, for every simple curve y c D between x and y in D, there exists a 
simple curve y' c D between x and y such that the rf//(y,}'') < 2 \fd/n and 
len c (y') < len c (y) + c/n for some c = c(d,D) > 0 (see Figure |4|. Since r] o u 
is continuous, y and y' are simple, and their Hausdorff distance is o(l), l |8.6| 
follows. 


Step II: Each geodesic (w.r.t. either d® ou or dist°) intersects only 0(n) cubes (in 
each of which the distortion is small). Consider a vertex of one of the cubes, 
and a facet of the same cube that does not intersect it. Since distances with 
respect to d, J0 „ are bounded from below by r] = min^ t](n) times the Euclidean 
distances, the distance between the vertex and the facet is at least rj/n. Therefore, 
a ball in (D, d^ ou ) of radius i//2/i intersects at most 2 d cubes. It follows that all 

geodesics in (D, d^ ou ) (or more generally, every curve of length diam c D or less) 
intersect at most 

diam c D d 2 d+l diam c D 

- 7 7T— ' 2 =-— • n 

Tj/2n i] 

of the cubes n H/ i ,□ n ,k n +m n (see Lemma 5.3 in lKM15bl for a similar argument). 
By the same reasoning, each geodesic in either (D, d^ ou ) or (D, dy ltl ) intersects at 
most 2 d+1 n diam e D/i] cubes. 

Next, cover similarly Dr by "cubes" luN = 7TR(n H<! ); also, define D R and Dr in 
an analogous way. Applying Proposition |8. 1 1 to each cube we obtain that 

lim P, (/ R(dis7TR| n „, < 3 VdA„/tt Vi = 1 ,.. . ,k n + m„) = 1 . (8.7) 

R —»oo 

For a given R , denote 

A r = {dis n R \a ni < 3 VdA„/« Vi = 1,... ,k„ + m„}. 


53 







Figure 4: Illustration of the construction of step I of the proof. The boundary of 
the domain D is marked in black, and the dotted lines represent the grid that 
forms The boundaries of D and D are marked with blue and red lines. 


respectively. For every simple curve y c D between two points in D there is a 
curve y' c D of similar length which is close to y in the Flausdorff distance (see 
( |8.6) ). Therefore, the inclusion of D in D has a vanishing distortion as n —» oo, 
as stated in <18.5 1. 
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Assume that Ar holds, and consider a vertex in one of the cubes in Dr, and a 
facet of the same cube that does not intersect the vertex. The distance between 
the vertex and the facet is at least (77 — 3 VrfA„)/n, which is positive for large 

enough n. By the same reasoning as above, each geodesic in either D ( , or Dr 
intersects at most 2 d+1 n diamD/(?] - 3 VdA„) cubes. 


Step III: Bounding the distortion. We now want to show that when Ar 
holds, i.e., the distortion within each cube is small, then the distortion of 
n r : (D,d^ ou ) —> (Dr, dist°) is small as well (and similarly with D). The idea, 
which is similar to the proof of Theorem 3.1 in I KM15 I. is that a geodesic y in 

(D, d® ou ) crosses 0(n) cubes at points x m , and at each crossing it accumu¬ 
lates a distance d 1]ou (xj,Xi + 1 ). When Ar holds, 

\d n o U (xi,x i+ {) - dist R (n R (x;),n R (x/+i))| < 3 ^ A ” , 

i.e. the distortion accumulated in each cube is 0(A n /n), so the total distortion is 
0(n) ■ 0(A n /ri) = 0(A„) —> 0. Here and below, the constants in O(-) only depend 
on d, diam t D and /], which are fixed throughout the proof. 

Formally, let x,y e D, and let y be a geodesic in (D, d® 0u ) between x and y. If 
y does not intersect the facets of the boxes transversely, we can take y to be 
infinitesimally longer such that it does. Therefore, without loss of generality 
y intersects the facets of the boxes at a finite number of points X \,... ,x m , with 
m < 2 d+1 n diam D/i] = 0(n). Denote x = xo and y = x m+ \. When the event Ar 
holds, 

dist 

< dqou( x i, x M) + m ■ 0(A„/«) = len d oJy) + 0(A„) 

i =1 n °“ 

= d ? ]au (x,y) + 0(A n ). 

A similar argument, on a geodesic a between 7 Xr(x) and nn(y), shows that 

tfou( x >y) ^ dist R (n R (x), n R (i/)) + 0(A„)- 
Therefore, we obtain that when Ar holds, 

sup | d° ou (x, y) - dist R (n K (x), n R (y))\ < CA„, ( 8 . 8 ) 

x,yeD 

for some C > 0 independent of n. Similarly, when Ar holds, 

sup | dyy oll (x, y) - dist^( 7 i R (x), tt r (i/))| < CA„. (8.9) 

x.yeD 


(n R (x),n R (y)) < dist R (7i R (xO,7i R (Xj + i)) < d t]OU (xi,x M ) 
i=l i= 1 V 


3 VdA„ 
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If D = D, then the proof is complete. The rest of the proof deals with the case 
where DcDcD. InD we have that d ROU < d ROU < dy OH , hence 

\d% u - <J < K u - d° oll \. (8.10) 

Similarly, distj^ < dist° < distjy in D R , hence 

| distjy - dist R | < | distj^ - dist^ |. (8.11) 

Using ( |8.5| , Q-fm), we obtain (similarly as in ( |8.4) ) that when A R holds, 
sup \d° ou (x,y) - dist R (7T R (x),7t R (l/))| 

x.yeD 

< sup | d° ou (x,y) - d° ou (x,y)\ + \d° oll (x,y) - distf(n K (x),n R (y))| 

x,yeD 

+ | distjy(n K (x), n R (y)) - dist°(n R (x),n R (i/))| 

< sup | d° ou (x,y) - df ou (x,y)\ + \d^ ou (x, y) - distf(7i R (x),7T R (y))| 

x,yeD 

+ | distf(n K (x), n R (y)) - dist^(7t R (x), tt r (t/))| 

< sup 2\d° ou (x,y) - d^ ou {x,y)\ + 2\ d^ ou (x,y) - distjy(7i R (x),n R (i/))| 

x,yeD 

+ \tfoJX'V) - dist^(n R (x), n R (y))\ 

< sup 2| d ROU (x,y) - d° ou (x,y)\ + 3CA„ = o(l), 

x.yeD 


where o(l) here is with respect to n and independent of R. The first inequality 
is a triangle inequality; the second follows from ( |8.10| ) and 1 8 . 11 ) ; the third is 
again a triangle inequality; and the last one follows from 18.81 and l|8.9). 


Since D is a 1/n-net in D and D R is a l/«-net in D R , we obtain that for tc r : 
(D,d^ ou ) —> (D r , dist R ), 

dis7i R = o(l). (8.12) 


Finally, let e > 0, and choose n large enough such that dis n R < e when A R 
holds. Note that this choice of n is independent of R. It follows that as R —» oo, 

P uA d GH((D R , dist R ),(D,„)) > e) < P„, R (dis n R > e) < P„ /R (Ay 0. 


8.2 Measure convergence 


In this section we prove part 2 of Theorem 2.2 
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For every e > 0, 


lim P„, R \3f e W(D), s.t. 

R—> oo \ 


fdy 

OOU 

Jd 




= o. 


where 

W(D) = {fe C(D) : ||/||„ < 1 and Lip(f) < 1}, 
cr(n) = e~ UKd/d , and kj is the volume of the d-dimensional unit ball. 

We start with a proposition which is a measure analog of Proposition |8.1| 

‘Proposition 8.2 Let D be a cube in ]R d and let u m j n = mino u > 0 and » max = 
maxo u < u,. Then 

hm P u,r (\v r (D r ) - p a0 u(D)\ > 5AcrLeb rf (D)) = 0, 
where Act = o(u min ) d - cr(u max ) rf - 


Proof: As in the proof of Proposition |8.1| consider the coupling between co u , 
and Since “min < u < u max , it follows that V R (",COu m J > v R (-; co u ) > 

v r(-; w >w)- Th en, 

\vr{Dr'/ (Vu) ~ Paoi/(Pl)| 

< | Vr(D r ;co u ) - v R (D R ;w Um J\ + \v R (D R ; a) Um J - p a ( Umax )(D)| + |/i CT ( Umax )(D) - p aou (D)\ 

< \v R (D R -,a) UmiD ) - v R (D R ; cu„ max )| + |vR(D K ;ai Umax ) - p a ( Umax ){D)\ + |pff( Umax )(D) - p a ( Umin )(D)\ 
— |vR(DR;ai Umln ) - p<j(ii mln )(D)| + 2 \v R (D R ;co Umax ) — Po(ti max )(D)| + 2 \pa(u mm )(D) — 

< | v R (D R ;w Ua J - p<j(u mln )(C>)| + 2 |v R (DR;a; limax ) - p a(Umax) (D)| + 2 Act Leb d (D), 

from which the result follows by Proposition |h4] ■ 

We now prove part 2 of Theorem |2.2| In order to simplify the notation, we will 
consider v R as a measure on D (assigning zero measure to S R ), rather than on 
D r . In this way there is no composition with tt^ 1 . 

As in the proof of the first part, partition D into cubes D„ r i ,and 
let D = [J k ff m “ We extend both p aou and v R to D by setting D \ D to be a 
null-set. 

Let e > 0 and choose n > Vd/e large enough such t hat A ct < e in each n D 
and Lebrf(D) < 2 Lebrf(D). By applying Proposition |8.2| to each we obtain 
that 

lim P„,r(Gr) = 1, 

R—*oo 

where 

Gr = {|vr(d„,;) — p 5e Leb^(n a ; ) Vi — 1, • • • ,k n + m n |. 
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Let / e W(D). Due to the choice of n and since Lip(/) < 1, it follows that 
sup n nD / - inf D)i jn o / < £ for every i. Denote f = inf D|ijn o /. Assume that Gr 
holds. Then, 


/"» k n ~\~TH n /~* 

I f d^ooii If dvR ^ ^ t I f djigou ~ I f dvR 

JD JD Jn nri Jn nii 


kn+m n 


+ \fi\ I flo ou 


r* 

— I f d^Qo u ~ fi liaou(j-^n,i) 

i =1 ^ a n,i 
k n +m n 

— L6b^(n W/Z ) + ||y||oo | l^oou(j^n,i) ~~ ^r(|—U, z)| 

z=1 

< 7c ^ Lebd(n„,0 = 7cLeb d (D) < 14cLeb d (D). 


fi vr(d 


n,i) - f 
dn„,j 


fdv R 


i=1 


Therefore, as R —* oo, 
Ph,r p/ e W(D), s.t. 


I /" dUijou I 
JD JD 


fdv R 


> 14c Leb d (D) < P U/R (G^) —» 0, 


which completes the proof. 


8.3 Asymptotic surjectivity 


In this section we prove Part 4 of Theorem 2.2 
For every c > 0, 


lim P„, R (d H (P, n-^D',,)) > c) = 0, (8.13) 

R —»oo ' ' 


where D R = {x e Dr : |tt r 1 (x)| = 1} and dn is the Hansdorjf distance in 

]R d . 


This is an immediate consequence of Proposition |7.2| using the above coupling 
betwe 
Then 


between nand u max = max D u. Denote D R = [x e D R (a> Um J : \n^(x-,oj Uia J\ = 1}. 


n-\D R ) c n R (D R ) c D. 

By Proposition j7?2] for every c > 0 

lim P R (d H (D,n^(D' R )) > c) = 0, 

R—>oo ' ' 

and combining ( |8.15) > with 18.141 we obtain | 8.13j . 


(8.14) 


(8.15) 
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8.4 Proof of Corollary |2.3| 

We follow here definitions and notations of Definition 12.II 

Since 7 Tr : D —» Dr is always defined on the whole space D and always onto 
Dr, Part 1 of Theorem |2.2| as stated in ( |8.1 ) immediately implies (JzTS}, which 
proves the first part of Corollary |2.3| 

For the second half of the first part, we need to prove that 


(i) n R ] : D', —> D is an e„-approximation 

(ii) For every / e C(D), 


f°< dv ^ = l/ d ^ 


(8.16) 


Note that since Vr(Dr \ D' R ) = 0 by definition, we can replace the integrals 
on the left-hand side by integrals on Dr. 

Part 1 (the event in | 2.15) ) implies that dis ti^ 1 < e„. Part 4 (the event in ( 2.17) ) 
implies that the e,,-neighborhood of n^(D' R ) with respect to the Euclidean 

metric in R rf contains D, hence also with respect to d ROU . By applying Hr„, 
it follows that the £ n -neighb orh ood of D', in D R n = 7iR n (D) with respect to 
is indeed the whole D 


R„ 


R n . This shows that tt^ 1 : D' R 


D is indeed an 


^-approximation. 

To prove ( 8.16) , note that it suffices to prove it for Lipschitz functions on D. 
Indeed, suppose we proved ( |8. 16) for Lipschitz functions, and let / e C(D). For 
every £ > 0 there is a Lipschitz function g on D such that ||/ - ylU < £, hence 


lim sup 


I f on R„ dv Rn- I fdy-a. 

Jd 

I g° n Rl dv R„- I S d l l 

Jd 


>Dr„ 

< lim sup 

n 

< 2e Lebrf(D). 


+ £Vr(DrJ + £p 

OOU (D) 


To prove ( 8.16) for a Lipschitz function /, note that for M > 0 large enough, 
f/M e W(D). Therefore, when Part 2 (the event in ( 2.16) ) holds. 


I f°^ R l dv R„- I f d l l a 

Jd r „ Jd 


< Me n —> 0, 


which completes the proof. 
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